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INTRODUCTION 


The  ultrasonic  vibration  potential,  originally  predicted  by  Debye1  in  1933,  refers  to 
voltage  generated  when  a  colloidal  or  ionic  solution  irradiated  with  ultrasound.  In  a  colloid,  the 
voltage  production  is  the  result  of  differential  motion  between  the  particles  in  suspension  and  the 
surrounding  fluid.  Colloids  are  suspensions  of  charged  particles  in  solution.  Some  examples 
include  milk,  muddy  water,  and  India  ink.  Normally  the  counter  charge  around  each  particle  is 
spherical;  however  when  a  sound  wave  propagates  through  the  suspension,  there  is  a  differential 
motion  of  the  particles  and  the  fluid  owing  to  their  different  inertias.  Typically,  the  particles  are 
dense,  and  tend  to  remain  stationary  in  space.  The  result  of  the  differential  motion  is  that  the 
fluid  oscillates  back  and  forth  over  the  surface  of  the  particles.  Since  the  counter  charge  moves 
with  the  fluid,  an  electric  dipole  is  created  at  the  site  of  each  particle  owing  to  the  distortion  of 
the  normally  spherical  charge  distribution.  Over  a  macroscopic  length,  the  polarization  from 
addition  of  the  dipoles  results  in  a  voltage  that  can  be  measured  with  a  pair  of  electrodes  in 
solution. 

We  have  found  that  whole  blood  generates  a  signal  that  is  substantial  compared  with 
animal  tissue.  The  potential  benefit  of  the  method  relative  to  conventional  x-ray  and  ultrasonic 
imaging  is  the  contrast  that  should  be  obtainable  in  a  vibration  potential  image  would  come  from 
the  presence  of  blood,  which  is  a  marker  for  many  invasive  tumors.  The  resolution  of  the 
method  is  governed  ultimately  by  the  wavelength  of  the  ultrasound.  In  the  experiments  carried 
out  with  one-dimensional  objects,  the  resolution  using  a  frequency  domain  version  of  the 
instrument,  resolution  much  higher  than  that  dictated  by  the  radiation  wavelength  is  possible. 

We  have  also  investigated  phase  contrast  x-ray  imaging  for  tumor  detection.  The 
principle  on  which  this  kind  of  imaging  relies  is  on  a  recording  of  the  phase  changes  that  x-rays 
experience  on  traversing  a  body.  In  principle,  phase  contrast  measurement  is  on  the  order  of 
1000  times  as  sensitive  as  absorption  contrast  for  most  of  the  useful  x-ray  region  of  the  spectrum. 
The  idea  of  the  method  is  old,  but  its  application  to  the  x-ray  region  of  the  spectrum  in  a 
convenient  geometry  was  first  reported  by  Wilkins  and  coworkers9.  To  observe  phase  contrast  in 
an  image  using  the  “in-line”  geometry  introduced  by  Wilkins,  the  x-ray  source  must  have  high 
spatial  coherence,  but  not  necessarily  high  temporal  coherence,  that  is,  there  is  no  severe 
consequence  to  using  polychromatic  x-rays.  X-ray  sources  that  have  been  developed  for  phase 
contrast  imaging  include  the  microfocus  x-ray  tube  and  the  synchrotron,  the  former  generating  an 
approximation  to  a  spherical  wave,  and  the  latter  a  plane  wave. 

In  the  in-line  geometry,  objects  with  differing  densities  induce  different  phase  lags  in  the 
x-rays  as  they  traverse  a  body  that  are  manifested  by  deflection  of  the  x-rays  and  the  formation 
of  interference  fringes  in  the  image.  The  peculiar  feature  of  phase  contrast  imaging  is  the 
highlighting  of  objects  with  different  densities  with  a  sharp  outline  in  the  image.  Our  method 
uses  ultrasonic  radiation  pressure  to  produce  displacements  in  an  object.  The  method  we  use  for 
generating  an  image  is  to  take  a  single  x-ray  image  using  a  CCD  camera  and  to  transfer  the 
image  to  a  computer.  Next  a  second  x-ray  image  is  made  but  in  the  presence  of  the  ultrasound. 
The  two  images  are  subtracted  pixel  by  pixel  in  the  computer  to  form  the  final  image.  The 
subtraction  results  in  a  heightening  of  the  phase  contrast  features  of  the  image  and  a  removal  of 
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the  absorption  features.  The  utility  of  the  method  for  tumor  imaging  is  that  the  tumor,  owing  to 
its  mechanical  properties  moves  as  a  unit  in  response  to  the  ultrasound.  Often  physicians  use  the 
same  mechanical  property  of  tumors  in  what  is  referred  to  as  “palpation”  to  locate  the  presence 
of  tumors  in  soft  tissue  through  use  of  finger  pressure.  The  perception  of  a  tissue  mass  that 
moves  as  a  unit  indicates  the  possible  presence  of  a  tumor.  Since  tumors  have  different  elastic 
properties  from  healthy  tissue,  their  sound  speed  differs  from  that  of  the  surrounding  tissue  so 
that  a  differential  radiation  force  is  experienced  by  the  tumor.  The  objective  of  the  approach 
under  investigation  here  is  to  combine  the  high  spatial  resolution  of  x-ray  mammography  with 
the  specificity  of  palpation  to  produce  an  image  that  highlights  objects  that  move  as  a  unit. 


BODY 

The  qualitative  explanation  given  above  for  generation  of  the  ultrasonic  vibration 
potential  by  a  colloid  has  been  investigated  by  a  number  of  authors  since  its  prediction  by  Debye. 
Reviews  have  been  published  by  Povey3,  Zana  and  Yeager4,  Babchin,  Chow  and  Sawatzky,5  and 
O’Brien,  Cannon,  and  Rowlands6'9.  The  theory  of  the  vibration  potential  for  dilute  suspensions 
of  O’Brien  and  coworkers6"9  is  more  rigorous  than  any  other  formulation,  and  is  widely  accepted 
by  workers  in  the  field.  According  to  O’Brien’s  theory,  the  electric  field  across  a  region  of  the 
colloid  generated  by  a  pressure  gradient  Vp  is  given  by, 

E  =  fAPm  ^  Vp,  (1) 

Pm  ^ 

where  A pm  is  the  difference  between  the  densities  of  particle  and  the  fluid,  f  is  the  volume 

fraction  of  the  particles,  pF  is  dynamic  mobility  of  the  particle,  cr*  is  the  conductivity,  and  pm  is 

the  fluid  density.  Vibration  potential  imaging  is  based  on  generation  of  a  voltage  across  a  pair  of 
electrodes  by  an  ultrasonic  wave  that  traverses  the  body.  The  important  part  of  the  research 
conducted  here  is  to  determine  the  presence  of  a  colloidal  or  ionic  fluid  within  a  body  and  to 
determine  its  distance  from  the  electrodes  and  its  geometry.  This  requires  a  complete  theory  of 
image  formation,  the  starting  point  for  which  is  Eq.  1 . 

The  work  with  phase  contrast  imaging  has  focused  on  setting  up  a  laboratory  provided 
with  a  microfocus  x-ray  tube,  a  liquid  nitrogen  cooled  CCD  camera,  fiber  optic  scintillator  plate, 
optics,  and  sample  mounts.  The  apparatus  is  interfaced  with  a  computer  so  that  the  x-ray  tube 
can  be  turned  on  and  off  remotely,  and  the  CCD  image  transferred  to  the  computer  for  storage. 
The  method  of  producing  an  image  in  these  experiments  is  to  take  a  single  phase  contrast  image 
of  the  body,  and  to  store  the  image  in  a  computer.  Second,  ultrasound  is  applied  to  the  body  that 
moves  objects  having  densities  or  sound  speeds  different  from  that  of  the  body.  The  resulting 
radiation  force  causes  a  movement  of  the  object  so  that  a  second  image  records  a  modified 
image.  The  two  x-ray  images  are  then  subtracted,  resulting  in  an  image  that  has  most  of  the 
absorption  component  of  the  image  removed,  but  with  a  highlighted  phase  contrast  component 
left  in  the  image. 

KEY  RESEARCH  ACCOMPLISHMENTS 

The  theory  of  production  of  a  vibration  potential  image  has  been  formulated  (see 
“Imaging  with  the  Ultrasonic  Vibration  Potential”:  A  Theory  for  Current  Generation”, 
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Ultrasound  in  Med  &  Biol.,  Appendix  I)  for  an  experimental  arrangement  where  parallel 
electrodes  are  placed  around  the  body,  and  the  sound  beam  propagates  perpendicularly  to  the 
plates.  The  object  is  considered  to  lie  within  a  limited  region  of  the  body,  and  the  dielectric 
constant  and  the  conductivity  of  the  body  are  considered  to  be  uniform.  The  main  conclusion  of 
the  theory  thus  far  formulated  is  that  the  current  in  the  electrodes  is  given  by 

I  _  i(p_  f  a(x,co)-@-p(x)dxdydz. 

h  J  v  oz  (2) 

which  can  be  expressed  as 


1  = 


-M.  1  \ya(x,co)]zp(x)dxdydz. 
h  J  v 


(3) 


where  a  is  the  distribution  of  the  colloid  in  space,  co  is  the  frequency  of  the  ultrasound,  h  is  the 
distance  between  the  electrodes,  and p  is  the  pressure. 


Equations  2  and  3  form  the  basis  for  vibration  potential  imaging.  For  the  first 
experiments  we  have  looked  at  straightforward  rastering  of  the  transducer  to  produce  amplitude 
images.  The  apparatus  is  shown  in  Fig  1  of  Appendix  II  (“Imaging  based  on  the  Ultrasonic 
Vibration  Potential”  Appl.  Phys.  Lett )  The  paper  describes  two  methods  of  imaging,  one  based 
on  time  of  appearance  of  the  vibration  potential,  and  a  second  method  based  on  changing  small 
changes  in  frequency.  Our  approach  thus  far  in  this  research  has  been  to  take  colloidal  objects 
that  have  specific  geometries  and  to  predict  the  frequency  dependence  of  the  in  phase  and 
quadrature  components  of  the  current  to  confirm  the  theory  and  to  explore  fundamental 
properties  of  the  imaging.  Figs  1  and  2  show  some  recent  images  taken  with  silica  phantoms  in 
agarose. 


Fig.  1  (Feft)  Photograph  of  a  top  view  of  a  two  silica  colloid  discs,  one  beneath  another.  A  side 
view  of  the  phantom  is  shown  in  Fig  2.  The  image  is  obtained  by  scanning  the  focused 
transducer  over  a  rectangle  30  mm  by  30  mm.  in  0.5  mm  steps.  The  signal,  after  amplification  in 
a  low  noise  preamplifier,  is  sent  to  a  lock-in  amplifier.  (Right)  Vibration  potential  image  of  the 
phantom. 
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Fig:  2  (Left)  Side  view  of  the  two  silica  colloid  disk  phantom.  The  2D  scanning  in  Fig.  1  can 
discern  the  top  disk  disc,  but  with  the  bottom  one  missing.  Determination  of  the  presence  of  the 
second  disk  is  carried  out  by  fixing  the  transducer  at  the  center  position  and  sweeping  the 
frequency.  By  fitting  the  in-phase  component  of  the  signal,  both  the  length  and  vertical  position 
of  each  disk  can  be  determined. 

We  have  proven  the  feasibility  of  our  method  of  subtractive  imaging  based  on  acoustic 
radiation  pressure  to  modify  phase  contrast  images.  See  Appendices  III  and  IV.  To  date,  our 
work  has  used  a  microfocus  x-ray  tube  and  a  CCD  camera  with  a  phosphor  screen  to  make  x-ray 
images  with  phase  contrast.  It  is  well-known  that  the  method  requires  spatial  coherence  in  the  x- 
ray  source,  as  is  produced  by  a  synchrotron  or  a  point  radiation  source.  We  are  developing  two 
approaches  for  providing  point  sources  of  radiation.  We  use  a  commercial  microfocus  x-ray  tube 
for  taking  x-ray  images  at  present,  and  have  been  developing  a  second  source  based  on  focusing 
a  fs  laser  to  generate  x-rays. 
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Figure  3:  (a)  Conventional  phase  contrast  image  of  a  3.05  mm  diameter  Teflon  sphere  suspended 
by  a  glass  fiber  in  water  showing  both  absorption  and  phase  contrast,  (b)-(d),  subtracted  x-ray- 
acoustic  phase  contrast  images  taken  with  10  kPa,  0.6  MPa,  and  1.1  MPa  average  acoustic  power 
applied  to  the  transducer. 


Results  of  using  the  method  on  a  Teflon  sphere  in  water  are  shown  in  Fig.  The  figures 
where  the  radiation  force  is  0.6  MPa  and  1.1  MPa  show  the  movement  of  the  sphere  by  the 
radiation  forces,  and  the  virtual  elimination  of  the  absorptive  component  of  the  x-ray  image. 


KEY  RESEARCH  ACCOMPLISHMENTS 

A  rigorous  theory  of  the  formation  of  the  vibration  potential  has  been  formulated.  The 
theory  has  been  applied  to  forming  images  of  several  objects.  The  results  of  this  theory  can  be 
applied  to  model  simple  objects  such  as  a  spherical  body,  a  cylindrical  body,  and  a  layer,  or,  in 
general,  to  any  arbitrary  object  within  a  body.  The  method  has  been  applied  to  form  images  of 
colloidal  objects  through  use  of  the  time  of  appearance  of  the  vibration  potential  signal  as  well  as 
a  simple  rastering  over  the  surface  of  an  object  to  give  an  “amplitude  map”  of  the  colloid  in 
space.  The  results  are  summarized  in  Appendices  I  and  II. 

Acoustically  modulated  x-ray  phase  contrast  images  have  been  taken  of  Teflon  spheres  in 
agarose  demonstrating  the  method,  (as  shown  in  Appendices  III  and  IV)  and  giving  a  general 
estimate  of  the  acoustic  radiation  pressure  required  to  produce  tissue  movement  of  sufficient 
amplitude  to  be  recorded  within  the  resolution  of  the  x-ray  tube  and  CCD  camera  used  here. 
Images  of  density  modified  agarose  (agarose  containing  colloidal  carbon)  in  chicken  flesh  have 
been  recorded.  See  Appendices  III  and  IV  for  details  of  the  experiments. 

REPORTABLE  OUTCOMES 

The  following  have  been  published  or  are  under  review: 

“Vibration  Potential  Imaging:  Theory  and  Preliminary  Results”,  with  A.  C.  Beveridge,  S.  Wang, 
and  V.  Gusev,  Proc  SPIE  2004;  5320:  95-100 

“Vibration  Potential  Imaging”  with  V.  Gusev,  A.  C.  Beveridge,  and  S.  Wang,  Proc.  NIST  Conf. 
on  Thermophysical  Properties,  Boulder,  2004 

“Imaging  based  on  the  Ultrasonic  Vibration  Potential”  with  A.  C.  Beveridge  and  S.  Wang,  ( Appl . 
Phys.  Lett.  85,  5466  (2004) 

“Acoustic  Radiation  Pressure:  a  “Phase  Contrast”  Agent  for  X-ray  Phase  Contrast  Imaging,  with 
C.  Bailat,  T.  Hamilton,  C.  Rose-Petruck  Appl.  Phys.  Lett  85  4517  (2004) 

“Imaging  with  the  Ultrasonic  Vibration  Potential:  A  Theory  for  Current  Generation”  with  V.  E. 
Gusev  Ultrasound  in  Med.  and  Biol.  31,  273  (2005) 
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“The  Theory  of  Ultrasonic  Vibration  Potential  Imaging”  A.  C.  Beveridge,  S.  Wang,  and  V.  E. 
Gusev,  Journal  de  Physique  IV,  accepted  for  publication) 

“Acoustically  Modulated  X-Ray  Phase  Contrast  Imaging”  with  T.  Hamilton,  and  C.  Bailat,  C. 
Rose-Petruck,  Phys.  Med.  Bio.  49,  4985  (2004) 

“Acoustically  Modulated  X-Ray  Phase  Contrast  and  Vibration  Potential  Imaging”  with  a.  C. 
Beverideg,  C.  J.  Bailat,  T.  J.  Hamilton,  S.  Wang,  C.  Rose-Petruck,  and  V.  E.  Gusev  (Proc.  SPIE 
2005,  in  press)  selected  as  the  best  conference  paper 


CONCLUSIONS 

An  operating  instrument  for  producing  images  using  the  ultrasonic  vibration  potential  has  been 
constructed.  The  theory  of  vibration  potential  imaging  has  been  given  for  an  arbitrary  object. 

The  method  has  been  applied  to  form  images.  Phase  contrast  images  using  a  microfocus  tube 
with  acoustic  modulation  of  objects  within  a  body  have  been  made.  Both  imaging  methods  show 
promise  for  tumor  detection. 
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Imaging  based  on  the  ultrasonic  vibration  potential 

Andrew  C.  Beveridge,  Shougang  Wang,  and  Gerald  J.  Diebolda) 

Department  of  Chemistry,  Brown  University,  Providence,  Rhode  Island  02912 
(Received  16  June  2004;  accepted  27  September  2004) 

An  ultrasonic  wave  traversing  a  colloidal  suspension  causes  distortion  of  the  charge  distributions  at 
the  sites  of  individual  colloidal  particles  producing  a  voltage  known  as  the  ultrasonic  vibration 
potential.  We  show  how  imaging  of  colloidal  regions  within  a  body  can  be  carried  out  using  a  beam 
of  ultrasound  to  produce  a  radio  frequency  vibration  potential.  A  theory  for  image  formation  shows 
that  Fourier  transformation  of  vibration  potential  signals  processed  by  a  mixer  and  low  pass  filter 
gives  the  spatial  distribution  of  colloid.  The  salient  feature  of  the  method,  insofar  as  medical 
imaging  is  concerned,  is  its  contrast  mechanism.  ©  2004  American  Institute  of  Physics. 
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Colloids  are  suspensions  of  charged  particles  in  a  liquid 
with  a  counter-charge  distributed  in  the  fluid  around  each 
particle.  The  counter-charge,  which  is  normally  a  spherical 
distribution  around  the  particles,  gives  the  solution  overall 
charge  neutrality  and,  together  with  the  charge  on  the  particle 
surface,  stabilizes  the  suspension  against  particle  agglomera¬ 
tion.  When  sound  propagates  through  a  suspension  where  the 
particles  have  either  a  higher  or  lower  density  than  that  of  the 
surrounding  fluid,  the  amplitude  and  phase  of  the  particle 
motion,  owing  to  the  difference  in  inertia  between  the  par¬ 
ticle  and  the  volume  of  fluid  it  displaces,  differs  from  that  of 
the  fluid  so  that  fluid  flows  back  and  forth  relative  to  the 
particle  on  alternating  phases  of  the  acoustic  cycle.  Since  the 
counter-charge  is  carried  by  the  fluid,  the  oscillatory  motion 
of  the  fluid  relative  to  the  particle  distorts  the  normally 
spherical  counter-charge  distribution  creating  an  oscillating 
dipole  at  the  site  of  each  particle,  which,  when  added  over  a 
half  wavelength  of  the  sound  wave,  results  in  a  macroscopic 
voltage  that  can  be  recorded  by  a  pair  of  electrodes  placed  in 
the  solution.  Similar  considerations  of  particle  inertia  show 
that  vibration  potentials  are  generated  in  ionic  solutions  as 
well.  The  theory  of  the  colloidal  vibration  potential1-3  gives 
the  electric  field  E  generated  by  an  acoustic  pressure  p  as 
E  =  -((pAppiE/ pK*)  Vp,  where  cp  is  the  volume  fraction  of 
particles,  Ap  is  the  density  difference  between  the  particles 
and  the  solvent,  ptE  is  the  electrophoretic  mobility  of  the 
particles,  p  is  the  solution  density,  and  K*  is  the  conductivity 
of  the  suspension.4  According  to  linear  acoustics  of  an  invis- 
cid  medium,  the  gradient  of  the  pressure  is  related  to  the 
fluid  acceleration  ii(r  ,t)  through  ii  =  -Vp/p,  which,  for  a 
plane  wave  propagating  along  the  z  axis,  can  be  written  in 
terms  of  the  displacement  u(z,t)  at  a  time  t  as  ii(z,f) 
=  -d)2u(z,t ),  where  o  is  the  radial  frequency  of  the  wave. 
Since  the  potential  across  the  region  of  integration  AT7'  is 
given  by  ATr(r)  =  /E-  dz,  it  follows  for  a  one-dimensional 
geometry5,6  that 

ATr(0  =  T  J  u(z,t)c(z)dz,  (1) 

where  the  quantities  in  the  integral  marked  with  a  hacek  are 
dimensionless,  i.e.,  z=z/L  and  u  =  u/u0 ,  u0  is  the  magnitude 
of  the  acoustic  displacement  in  the  plane  wave,  L  is  the 
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length  of  the  colloidal  region,  and  T  is  a  proportionality  con¬ 
stant  given  by  r  =  -(D2cpApLEu0L/ pK*,  with  the  integration 
extending  over  the  spatial  extent  of  the  sound  wave. 

The  apparatus  used  for  vibration  potential  imaging,  as 
shown  in  Fig.  1,  consisted  of  a  cell  containing  the  object  of 
interest,  equipped  with  acoustic  delay  lines  on  the  top  and 
bottom  of  the  cell,  with  two  electrodes  at  opposite  ends  of 
the  cell,  mounted  directly  on  the  faces  of  the  delay  lines.  One 
electrode  was  attached  to  an  aluminum  housing  surrounding 
the  cell,  which  forms  a  ground,  and  the  other  electrode  was 
connected  to  a  feedthrough  and  cable  leading  to  a  rf  ampli¬ 
fier  (Analog  Devices,  Inc.  Model  AD  8021)  with  a  voltage 
amplification  of  100.  Ultrasound  was  generated  by  a 
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FIG.  1.  Diagram  of  the  apparatus.  Ultrasound  from  the  transducer  enters  and 
exits  the  colloidal  region  of  the  cell  through  a  variable  length  (approxi¬ 
mately  50  mm  long)  acoustic  delay  lines.  An  aluminum  housing,  which 
serves  as  the  electrical  ground,  supports  the  cell  and  the  delay  lines.  The 
delay  lines  were  either  acrylic  blocks  with  copper  tape  electrodes  mounted 
on  the  faces  of  the  delay  lines  or  water  columns  with  90  gm  wires  used  as 
electrodes  and  plastic  membranes  used  to  separate  the  colloid  from  the 
water.  For  experiments  where  an  image  is  generated,  the  transducer  is 
scanned  in  the  XY  plane  and  the  vibration  potential  signal  after  processing  is 
stored  in  a  computer  as  a  function  of  position.  Inset:  Oscilloscope  trace  of 
vibration  potential  (in  arbitrary  units)  vs  time  generated  by  approximately 
six  pulses  of  1  MHz  ultrasound  traversing  a  water  delay  line  into  a 
2.538-cm-thick  layer  of  colloidal  silica. 
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FIG.  2.  Oscilloscope  traces  of  the  vibration  potential  generated  by  32  pulses 
of  (a)  1.0165  MHz,  (b)  1.0069  MHz,  and  (c)  0.9865  MHz  ultrasound  tra¬ 
versing  a  saline  solution  delay  line  into  2. 5 -cm- thick  layer  of  colloidal 
silica.  Signals  were  averaged  over  64  repetitions  of  the  pulse  burst.  The 
oscilloscope  voltage  is  referred  to  the  electrodes  after  amplification.  The 
initial  (55-75  juls)  and  final  (85-105  jus)  parts  of  the  wave  forms  corre¬ 
spond  to  entrance  of  the  ultrasound  into  the  colloidal  layer  and  exit  from  the 
layer,  respectively.  The  central  portions  of  the  wave  forms  show  the  effects 
of  integration  of  the  wave  over  the  entire  length  of  the  colloidal  layer.  For 
the  0.9865  MHz  trace,  the  colloidal  layer  thickness  is  an  integral  number  of 
sound  wavelengths  giving  complete  cancellation  of  the  voltage  in  the  center 
of  the  wave  form. 


2.54-cm-diameter  LiNb03  transducer  driven  by  a  program¬ 
mable  pulse  generator,  or  by  a  function  generator  that  drove 
a  power  amplifier. 

An  experiment  was  carried  out  where  a  short  burst  of 
ultrasound  consisting  of  six  cycles  of  a  1  MHz  sine  wave, 
generated  at  a  repetition  rate  of  60  Hz,  was  applied  to  the 
cell  which  contained  a  2.538-cm-thick  colloidal  layer.  The 
spatial  extent  of  the  pulse  burst  was  smaller  than  the  length 
of  the  colloid  layer  in  the  cell.  The  signal  from  a  colloidal 
silica  sample7  recorded  by  the  oscilloscope  is  shown  in  the 
inset  of  Fig.  1.  As  the  ultrasound  burst  either  enters  or  leaves 
the  colloid,  an  alternating  vibration  potential  is  produced. 
However,  when  the  entire  six  cycles  of  the  burst  of  ultra¬ 
sound  are  totally  confined  within  the  colloidal  region,  the 
positive  and  negative  contributions  to  the  potential,  in  agree¬ 
ment  with  Eq.  (1),  integrate  to  zero  giving  a  null  voltage. 

For  bursts  of  ultrasound  whose  extent  in  space  is  larger 
than  the  thickness  of  the  colloidal  layer,  the  voltages  from 
individual  half-cycles  of  the  sound  wave,  according  to  Eq. 
(1),  can  be  nonzero  or  add  to  zero  depending  on  thickness  of 
the  layer  and  the  wavelength  of  sound.  Figure  2  shows  vi¬ 
bration  potential  signals  at  three  different  frequencies  pro¬ 
duced  by  a  32  cycle  burst  of  ultrasound  whose  spatial  extent 
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is  large  compared  with  the  2.538  cm  layer  thickness,  show¬ 
ing  the  effect  of  voltage  integration  over  the  length  of  the 
cell.  Evaluation  of  Eq.  (1)  for  a  wave  that  fills  the  colloidal 
region  shows  that  the  zeroes  in  the  signal  occur  when  nX 
=L ,  where  n  is  an  integer,  X  is  the  wavelength  of  the  sound, 
and  L  is  the  length  of  the  layer.  It  follows  for  two  frequencies 
f2  and  fi,  which  have  zeroes  for  integer  values  of  the  wave¬ 
length  n2  and  nx,  that  (n2-ni)c=L(f2-fi)  where  c  is  the 
sound  speed  which  permits  determination  of  the  layer  thick¬ 
ness.  Data  taken  over  a  more  extensive  range  than  shown  in 
Fig.  2,  from  0.8984  to  1.0765  MHz,  gave  the  length  of  col¬ 
loid  as  2.553  +  0.025  cm. 

There  are  three  kinds  of  imaging  that  follow  from  Eq. 
(1)  for  determining  the  Z  dependence  (as  defined  in  Fig.  1)  of 
colloidal  distributions  in  space  for  a  focused  transducer 
scanned  in  the  XY  plane.  The  first  method  uses  the  amplitude 
of  the  vibration  potential  recorded  with  a  lock-in  amplifier  as 
the  transducer  is  rastered  in  the  XY  plane  to  form  an  image. 
The  second  method,  analogous  to  conventional  pulse-echo 
imaging,  is  based  on  recording  the  delay  time  between  the 
appearance  of  the  vibration  potential  signal  and  the  time  of 
launching  of  the  acoustic  burst.  The  delay  time,  together  with 
the  sound  speed,  gives  the  distance  between  the  colloidal 
object  and  the  transducer.  For  the  third  method,  consider  an 
arbitrary  distribution  of  colloid  in  space  c(z )  traversed  by  a 
plane  wave  propagating  in  the  Z  direction.  It  is  straightfor¬ 
ward  to  show  that  accumulation  of  in-phase  and  quadrature 
voltages  from  a  mixer  and  low  pass  filter,  recorded  over  a 
range  of  k ,  or,  equivalently,  over  a  range  of  frequencies, 
gives  the  distribution  of  colloid  in  space,  c(z),  through  Fou¬ 
rier  transformation  as 


4 

c(z)  =  - 


J  Lj(k)cos(kz)dk  +  J  Lg(k)sin(kz)dk 


(2) 


where  L/(k)  and  Lqik)  are  the  in-phase  and  quadrature  com¬ 
ponents  of  the  signal  following  low  pass  filtering.  Equation 
(2)  shows  that  the  vibration  potential  can  be  inverted  to  de¬ 
termine  any  arbitrary  colloid  distribution  in  space  c(z).  It  is 
noteworthy  that  a  layer  of  uniform  thickness  extending  along 
the  z  axis  a  distance  L  from  the  origin  gives  Lj(k) 
=  (r/2&)sin (kL)  and  Lg(k)  =  (r/2k)[l -cos(kL)],  the  first  of 
which  can  be  written  in  terms  of  the  familiar  sine  function  as 
L7(k)  =  (rL/2)sinc(^L).  The  amplitudes  of  the  central  por¬ 
tions  of  the  traces  in  Fig.  2  give  three  points  on  the  sine 
curve  Lj(k). 

Experiments  were  carried  out  with  objects  to  give  a 
demonstration  of  vibration  potential  imaging.  An  agarose 
block  containing  colloidal  discs  was  placed  at  the  bottom  of 
the  cell  with  the  agarose  making  contact  with  the  electrode 
leading  to  the  preamplifier.  The  apparatus  shown  in  Fig.  1 
was  modified  so  that  a  saline  solution  in  a  cylindrical  enclo¬ 
sure  replaced  the  acrylic  delay  line  at  the  top  of  the  cell.  The 
salt  solution,  which  formed  one  of  the  electrodes  for  the  cell, 
made  contact  with  the  metal  housing  of  the  cell  providing 
both  a  ground  for  the  rf  signal  and  a  delay  line  for  a  beam  of 
ultrasound8  from  a  5.1-cm-diam,  15.2  cm  focal  length 
LiNb03  transducer  (Valpey  Fisher  Model  El  129).  The  ultra¬ 
sound  beam  was  approximately  3  mm  in  diameter  beam  over 
the  length  of  the  sample.  The  transducer  was  driven  by  a 
pulse  generator  that  produced  900  kHz  voltage  bursts  28  /xs 
long  with  a  repetition  rate  of  50  Hz,  amplified  by  a  gated 

power  amplifier.  The  position  of  the  transducer  relative  to  the 
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FIG.  3.  (Color)  (a)  Photograph,  (b)  two-dimensional  amplitude,  and  (c) 
three-dimensional  surface  vibration  potential  images  of  thin  layer  colloidal 
phantoms  cast  in  agarose.  In  (c),  plotted  with  the  ultrasound  beam  propa¬ 
gating  upwards,  the  depth  was  determined  from  the  delay  time  of  the  vibra¬ 
tion  potential  signal  relative  to  the  launching  of  the  acoustic  burst;  only  the 
upper  surfaces  of  the  colloidal  objects  are  recorded.  The  XY  scan  data  were 
recorded  at  intervals  of  1  mm  with  an  integration  time  at  each  pixel  of  1  s. 
The  driving  voltage  on  the  focused  LiNb03  transducer  was  250  V 
peak-to-peak. 


phantoms  in  the  cell  was  scanned  by  precision  translation 
stages  controlled  by  a  computer.  Signals  from  the  preamp¬ 
lifier  were  fed  to  both  a  digitizer  and  a  lock-in  amplifier  as 
transducer  was  rastered  in  the  X  and  Y  directions  across  the 
sample. 

Figure  3  shows  a  photograph,  an  amplitude  image,  and  a 
time  of  arrival  image  of  a  25-mm-thick  agarose  block,  which 
contained  two  2-mm-thick  silica  colloid  discs,  13  mm  in  di¬ 
ameter,  made  by  mixing  silica  colloid  with  agarose,  with  one 
disc  located  11  mm  above  the  other  and  displaced  in  the  Y 
direction.  The  time  of  arrival  image  was  made  by  setting  a 
threshold  signal  to  yield  a  unit  amplitude  output,  and  dis¬ 
carding  all  signals  with  smaller  amplitudes  than  the  thresh¬ 
old.  Note  that  the  signal-to-noise  ratio  in  the  images  could  be 
improved  significantly  by  employing  an  ultrasonic  wave 
train  with  a  higher  duty  cycle  than  the  0.2%  duty  cycle  used 
here. 

Experiments  were  also  carried  out  to  determine  the  mag¬ 
nitudes  of  vibration  potentials  from  colloidal  and  ionic  solu¬ 
tions,  and  biological  samples.  The  apparatus  shown  in  Fig.  1 
was  equipped  with  acrylic  delay  lines,  and  a 
2.54-cm-diameter,  planar  pzt  transducer  driven  by  a  single, 
negative  voltage  spike  from  a  pulse  generator  (Panametrics, 
Inc.  Model  5058PR).  A  number  of  common  colloidal  and 
ionic  solutions9  were  shown  to  produce  sizeable  vibration 
potentials  (given  in  parentheses),  including  colloidal  gold 
(15  julY),  India  ink  (44  /zV),  colloidal  silver  (13  /zV),  LiCl 
(6.9  /zV),  NaCl  (23  /zV),  KC1  (39  /zV),  RbCl  (85  /zV),  and 
CsCl  (132  /zV).  Insofar  as  tissue  imaging  is  concerned, 
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muscle  tissue  from  chicken  breast,  beef,  and  pork  all  pro¬ 
duced  vibration  potentials  smaller  than  0.02  /zV  while  whole 
blood,10  which  is  both  colloidal,  as  a  result  of  the  presence  of 
red  blood  cells,  and  ionic,  from  dissolved  electrolytes,  gave 
comparatively  large  signals,  on  the  order  of  10  /zV  at 
500  kHz.  Further  experiments  with  colloids  placed  between 
4-mm-thick  layers  of  chicken  breast,  with  the  breast  tissue 
placed  in  contact  with  the  electrodes,  showed  that  vibration 
potential  generated  by  the  colloid  within  the  weakly  conduct¬ 
ing  muscle  tissue  is  transmitted  to  the  surface  electrodes  and 
can  be  recorded. 

While  the  resolution  in  both  vibration  potential  and  con¬ 
ventional  ultrasonic  pulse-echo  imaging  are  governed  by  the 
same  principles  of  linear  acoustics,  the  contrast  mechanisms 
through  which  the  images  are  formed  in  the  two  imaging 
methods  are  altogether  different.  In  vibration  potential  imag¬ 
ing,  contrast  is  determined  by  differences  in  electroacoustic 
properties  within  the  irradiated  body,  i.e.,  by  differences  in 
the  colloidal  or  ionic  concentration  or  even  the  presence  of 
such  solutions,  whereas  in  conventional  ultrasonic  pulse- 
echo  imaging,  contrast  is  determined  by  relative  differences 
in  acoustic  impedance.  In  comparing  vibration  potential  im¬ 
aging  with  established  methods  such  as  optical  or  x-ray  im¬ 
aging,  it  is  important  to  note  that  resolution  is  ultimately 
limited  by  the  wavelength  of  the  radiation  employed.11  An 
additional  factor  in  considering  the  resolution  of  vibration 
potential  imaging  is  that  reflections  of  the  acoustic  waves  at 
interfaces  can  cause  changes  in  the  vibration  potential  signal, 
complicating  the  analysis  of  data.  Insofar  as  tissue  imaging  is 
concerned,  the  large  difference  in  signal  amplitudes  between 
blood  and  muscle  tissue  found  here,  which  is  greater  than  a 
factor  of  500,  points  up  its  contrast  mechanism,  which  differs 
fundamentally  from  that  for  radiographic,  photoacoustic,  op¬ 
tical,  or  ultrasonic  imaging,  and  which  suggests  its  use, 
among  other  things,  for  detection  of  tumors,  which  are 
known  to  be  highly  vascularized. 

The  experiments  reported  here  were  supported  by  the  US 
Army  Medical  Research  and  Materiel  Command  under  Grant 
DAMD 17-02- 1-0307.  Opinions,  interpretations,  conclusions 
and  recommendations  are  those  of  the  author  and  are  not 
necessarily  endorsed  by  the  US  Army.  The  authors  are  grate¬ 
ful  for  the  assistance  of  S.  Woythayler  for  design  of  the  rf 
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*R.  Zana  and  E.  Yeager,  in  Modern  Aspects  of  Electrochemistry ,  edited  by 
J.  O.  Bockris,  B.  E.  Conway,  and  R.  E.  White  (Plenum,  New  York,  1982). 
2R.  j.  Hunter,  Colloids  Surf.,  A  141,  37  (1998). 

3A.  S.  Dukhin  and  P.  J.  Goetz,  Colloids  Surf.,  A  192,  267  (2001). 

4We  use  Eq.  (4.2)  from  R.  W.  O’Brien,  J.  Fluid  Mech.  190,  71  (1988). 
5For  a  more  complete  description  of  vibration  potential  imaging  in  three 
dimensions  see  V.  E.  Gusev  and  G.  J.  Diebold  (unpublished). 

6The  analysis  here  has  ignored  the  effects  of  reflection  of  ultrasound  at 
boundaries  within  the  irradiated  body. 

7The  colloidal  silica  Snow  Tex  (Nisson,  Inc.,  Houston,  TX)  was  an  aqueous 
suspension  of  20-30  nm  silica  with  a  particle  concentration  of  1015/cm3. 
8The  wave  front  in  a  focused  sound  beam  has  curvature  which  must  be 
taken  into  account  in  a  precise  formulation  of  the  problem. 

9The  colloidal  gold  (Ted  Pella,  Inc.,  Redding,  CA)  was  an  aqueous  suspen¬ 
sion  of  100-nm-diam  particles,  at  a  concentration  of  5  X  109/cm3.  India 
ink  (Higgins,  Inc.,  Levenburgh,  TN)  is  an  aqueous  suspension  of  carbon 
with  a  surfactant  added  for  stability.  The  colloidal  silver  was  an  aqueous 
suspension  of  80-nm-diam  particles,  at  a  concentration  of  1.1  X  1010/cm3. 
The  salt  concentrations  were  4  M.  Production  of  a  vibration  potential  by 
these  solutions  is  well  known  from  previous  reports  in  the  literature. 
10Whole  canine  blood  preserved  with  EDTA. 

nThe  resolution  of  nuclear  magnetic  resonance  imaging  is  not  determined 
by  the  wavelength  of  the  radiation  employed. 
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Acoustic  radiation  pressure:  A  “phase  contrast”  agent  for  x-ray  phase 
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We  show  that  the  radiation  pressure  exerted  by  a  beam  of  ultrasound  can  be  used  for  contrast 
enhancement  in  high-resolution  x-ray  imaging  of  tissue  and  soft  materials.  Interfacial  features  of 
objects  are  highlighted  as  a  result  of  both  the  displacement  introduced  by  the  ultrasound  and  the 
inherent  sensitivity  of  x-ray  phase  contrast  imaging  to  density  variations.  The  potential  of  the 
method  is  demonstrated  by  imaging  microscopic  tumor  phantoms  embedded  into  tissue  with  a 
thickness  typically  presented  in  mammography.  The  detection  limit  of  micrometer  size  masses 
exceeds  the  resolution  of  currently  available  mammography  imaging  systems.  The  directionality  of 
the  acoustic  radiation  force  and  its  localization  in  space  permits  the  imaging  of  ultrasound-selected 
tissue  volumes.  The  results  presented  here  suggest  that  the  method  may  permit  the  detection  of 
tumors  in  soft  tissue  in  their  early  stage  of  development.  ©  2004  American  Institute  of  Physics. 
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Phase  contrast  imaging  differs  from  conventional  x-ray 
shadowgraphy  in  the  mechanism  of  contrast  generation: 
while  conventional  shadowgraphy  depends  on  absorption  of 
x  rays,  phase  contrast  imaging  is  based  on  phase  changes  as 
x  rays  traverse  a  body  giving  rise  to  wave  interference  that 
results  in  intensity  changes  in  the  image.  Fresnel-Huygens 
theory,1’2  which  governs  image  formation  in  general,  de¬ 
scribes  contrast  in  an  image  as  dependent  on  both  absorption 
and  phase  contrast.  However,  only  when  sources  with  high 
spatial  coherence,  such  as  synchrotrons,3-5  microfocus  x-ray 
tubes,6’7  or  laser  plasma  x-ray  sources8’9  are  employed  is  the 
phase  contrast  component  of  the  image  visible.  Coherent 
sources  give  an  image  of  a  dense  object  embedded  in  lower 
density,  soft  tissue  as  a  shadow  of  the  object,  caused  by 
absorption,  surrounded  at  its  perimeter  by  light  and  dark  in¬ 
terference  fringes  arising  from  rapid  phase  variations  in  the 
radiation  at  the  interface  between  the  two  media.  A  detailed 
analysis  of  the  relative  contrast  produced  in  an  image  for  a 
fixed  difference  in  density  shows  that  phase  contrast  is  more 
sensitive  than  absorption  throughout  most  of  the  spectral  re¬ 
gion  commonly  used  for  diagnostic  x-ray  imaging.10  For  in¬ 
stance,  the  thickness  of  water  needed  to  produce  a  1  %  phase- 
contrast  at  36  keV  photon  energy  is  2500  times  smaller  than 
that  required  to  produce  the  same  contrast  though 
absorption.10 

Phase  contrast  in  an  image  can  be  markedly  enhanced  by 
recording  two  images  of  an  object,  one  where  the  object  is 
displaced  slightly  through  acoustic  radiation  pressure,11-13 
and  a  second  where  the  object  is  unmoved,  followed  by  sub¬ 
traction  of  the  two  images.  Acoustic  radiation  forces14-16  can 
arise  from  either  reflection  of  sound  by  an  object  or  as  a 
result  of  acoustic  impedance  changes  (i.e.,  a  variation  in  ei¬ 
ther  density  or  sound  speed)  or  by  dissipative  processes 
within  the  object  resulting  in  absorption  of  the  ultrasound 
energy.  A  number  of  “elastographic”  studies,  where  tissue 
movement  is  monitored  in  response  to  radiation  or  mechani¬ 
cal  pressure,11-13  have  shown  that  variations  in  Young’s 
modulus,  which,  in  part,  determines  sound  speed,  permits 
acoustic  differentiation  of  tissue.  It  is  known,  for  instance, 
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that  breast  tumors  have  a  significantly  different  Young’s 
modulus  from  surrounding  tissue,  and  that  tumors  move  as 
rigid  bodies  in  response  to  acoustic  radiation 
pressure.11-13'17’18 

The  subtractive  phase  contrast  imaging  method  de¬ 
scribed  here  is  based  on  acoustic  radiation  pressure  to  gen¬ 
erate  differential  movement  of  components  of  a  body  having 
variations  in  acoustic  impedance.  First,  a  reference  x-ray 
phase-contrast  image  is  taken  with  a  CCD  camera  and  stored 
in  a  computer;  next,  a  second  x-ray  image  is  recorded  while 
the  body  is  irradiated  with  a  directed  beam  of  ultrasound, 
under  the  same  imaging  conditions  as  the  reference  image. 
By  subtracting  the  two  images,  a  difference  image  of  the 
temporally  averaged  displacement  induced  by  the  ultrasound 
is  obtained. 

As  shown  in  Fig.  1,  x  rays  were  generated  by  a  micro¬ 
focus  x-ray  tube  (Oxford  Instruments,  model  UB-M1),  oper¬ 
ated  at  90  kV,  which  has  a  minimum  source  diameter  of 
25  jam  as  measured  using  a  knife  edge  imaging  procedure. 
The  x  rays  were  detected  by  a  Gd202S(Tb)  fiber  optic  scin¬ 
tillation  plate  (Hamamatsu,  Inc.,  model  J6676)  the  fluores¬ 
cence  from  which  was  imaged  onto  a  liquid  nitrogen  cooled 
CCD  camera  (Roper  Scientific  model  7382-0001)  interfaced 
to  a  computer.  The  source  to  detector  distance  was  2.6  m, 
while  the  sample  to  source  distance  was  35  cm,  providing  a 
magnification  of  7.5.  The  x-ray  beam  exiting  the  sample 
propagated  in  a  helium  atmosphere  in  order  to  reduce  x-ray 
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FIG.  1.  Diagram  of  the  experimental  apparatus.  X  radiation  generated  by  a 
microfocus  tube  penetrates  a  sample  and  is  detected  by  a  CCD  camera  that 
views  a  scintillator  plate.  The  sample  is  located  in  a  3-cm-diam  tube  filled 
with  water  coupled  to  an  ultrasound  transducer.  The  intensity  maximum  in 
the  x-ray  spectrum  detected  at  the  scintillation  screen  was  at  40  keV  owing 
to  x-ray  spectrum  hardening  in  the  sample. 
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FIG.  2.  (a)  Conventional  phase  contrast  image  of  a  3.05-mm-diam  Teflon 
sphere  suspended  by  a  glass  fiber  in  water  showing  both  absorption  and 
phase  contrast,  (b)-(d)  Subtracted  x-ray-acoustic  phase  contrast  images 
taken  with  10  kPa,  0.6  MPa,  and  1.1  MPa  average  acoustic  power  applied  to 
the  transducer. 


scatter.  The  sample  was  situated  in  the  middle  of  a  3 -cm 
-diam  PVC  tube  which  had  portions  machined  out  and  re¬ 
placed  with  Mylar.  One  end  of  the  sample  tube  was  termi¬ 
nated  by  a  synthetic,  sound-absorbing  material  while  the 
other  was  sealed  to  a  2-cm-diam,  1  MHz,  LiNb  transducer, 
located  3  cm  from  the  part  of  the  sample  traversed  by  the 
x-ray  beam.  The  transducer  was  driven  by  a  pulse  train  at  a 
repetition  rate  of  133  Hz  consisting  of  15-/zs-long  1  MHz 
bursts  amplified  by  a  power  amplifier  with  a  peak  power  of 
1  kW. 

The  principle  of  the  method  was  demonstrated  by  taking 
images  of  a  3  mm  diameter  Teflon™  bead  attached  to  a  glass 
fiber,  suspended  in  water,  and  irradiated  with  pulsed  ultra¬ 
sound.  A  conventional  phase  contrast  image  showing  both 
absorption  and  phase  contrast  at  the  perimeter  of  the  bead  is 
shown  in  Fig.  2(a).  Figures  2(b)-2(d)  show  subtracted  im¬ 
ages  taken  over  a  range  of  successively  higher  acoustic  pow¬ 
ers.  The  absence  of  features  in  Fig.  2(b)  shows  that  at  the 
lowest  acoustic  intensity,  the  displacement  of  the  bead  is 
smaller  than  the  resolution  of  the  x-ray  imaging  apparatus.  In 
Fig.  2(d)  the  displacement  of  the  bead19  with  an  acoustic 
pressure  of  1.1  MPa  was  determined  from  the  x-ray  image  to 
be  200  jum.  Two  effects  from  the  subtraction  of  the  image  of 
the  displaced  bead  from  one  of  a  stationary  bead  can  be  seen 
from  examination  of  Fig.  2(c)  and  2(d):  first,  the  absorptive 
component  of  the  contrast  has  been  removed  by  the  subtrac¬ 
tion  leaving  only  phase  contrast  in  the  image,  and,  second, 
the  vector  character  of  the  image  is  evident  in  that  interfer¬ 
ence  fringes  that  lie  parallel  to  the  direction  of  the  ultrasound 
induced  motion  are  eliminated  while  those  perpendicular  to 
the  direction  of  motion  are  highlighted. 

The  capability  of  the  method  for  distinguishing  only 
those  objects  irradiated  with  ultrasound  was  tested  by  imbed¬ 
ding  two  Teflon  beads  in  agarose  and  directing  a  focused 
beam  of  ultrasound  to  irradiate  only  the  bead  nearest  the 
x-ray  source.  Figures  3(a)  and  3(b)  show  a  conventional 
phase  contrast  image  and  an  acoustically  modulated  phase 
contrast  image,  respectively.  The  ultrasound  irradiated  bead 
is  easily  identified  with  virtually  no  trace  of  the  bead  located 
outside  the  ultrasound  field. 

In  order  to  evaluate  the  efficacy  of  applying  the  imaging 
method  to  tissue-like  samples,  the  displacements  of  micron¬ 
sized  pieces  of  agarose  embedded  in  chicken  breast  were 
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FIG.  3.  (a)  Top  view  of  the  sample  consisting  of  two  Teflon  beads  cast  in 
agarose.  The  two  beads  are  separated  by  3  cm.  (b)  Conventional  phase  con¬ 
trast  image  of  two  Teflon  beads  in  agarose  and  (c)  an  acoustically  modulated 
phase  contrast  image  showing  only  the  bead  irradiated  with  ultrasound.  The 
beads  were  located  2.5  cm  apart. 


measured.  The  agarose  was  enriched  with  carbon  nanopar¬ 
ticles  to  achieve  a  density  approximately  14%  higher  than 
that  of  the  chicken  breast  in  order  to  simulate  the  density 
difference  between  malignant  and  healthy  human  breast 
tissue.  ’  Figure  4  shows  a  series  of  images  of  the  chicken 
breast  with  the  embedded  objects.  The  flesh  included  fat 
globules  which  were  deflected  by  the  ultrasound  field  along 
with  the  agarose  objects.  All  images  contain  identical  de¬ 
tected  x-ray  doses  and  are  reproduced  with  identical  grey 
scale  ranges  of  ±5%  of  the  average  image  intensity.  The 
images  and  inserts  display  the  same  sample  regions  around 
three  selected  agarose  objects,  labeled  1-3  with  estimated 


sizes  of  40X53  jum  (1),  30X46  jum  (2),  and  20X26  /um 
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FIG.  4.  X-ray  images  of  carbon-enriched  agarose  objects  embedded  in 
chicken  flesh:  The  round  objects  are  fat  globules  in  the  chicken  flesh.  The 
entire  sample  is  immersed  in  water,  (a)  X-ray  image  taken  with  an  enlarged, 
75-/um-diam  x-ray  source;  (b)  conventional  phase  contrast  image  taken  with 
a  25-/xm-diam  x-ray  source,  and  (c)  x-ray  ultrasound  enhanced  x-ray  phase 
contrast  image  taken  with  25-/xm-diam  x-ray  source.  Tissue  regions  that  are 
not  shifted  by  the  ultrasound  field  are  indicated  by  arrows.  Inset:  magnified 
images  of  objects  1-3.  The  background  noise  remaining  after  subtraction  in 
image  (c)  is  the  shot  noise. 

(3).  Figure  4(a)  shows  an  image  of  the  chicken  breast  where 
the  x-ray  spot  size  in  the  microfocus  tube  was  deliberately 
enlarged  to  75  jum.  The  phase  contrast  features  of  the  image 
are  almost  entirely  eliminated  by  convolution  of  the  fringes 
over  the  spot  size  of  the  x-ray  source — only  the  absorption 
component  of  the  contrast  remains  in  the  image.  Fatty  mate¬ 
rial,  which  appears  as  nearly  circular  objects  in  the  image,  is 
visible  but  not  well  defined.  Three  different  agarose  pieces, 
referred  to  here  as  “objects,”  are  in  the  field  of  view  of  the 
imaging  system.  While  object  1  is  visible,  object  2  is  poorly 
defined,  and  object  3  is  essentially  indistinguishable  from  the 
background.  The  overall  visibility  of  the  objects  in  the  tissue 
is  low.  Figure  4(b)  shows  a  conventional  phase  contrast  im¬ 
age  with  the  x-ray  tube  operating  with  its  minimum  source 
diameter,  25  jum.  The  image  shows  the  fatty  masses  as  being 
clearly  outlined,  permitting  visual  identification.  Both  ob¬ 
jects  1  and  2  are  visible,  but  object  3,  again,  is  not  easily 
distinguished  from  the  background.  The  ultrasound  contrast- 
enhanced  x-ray  image  is  shown  in  Fig.  4(c).  The  ultrasound 
intensity  was  adjusted  to  give  a  tissue  displacement  of  ap¬ 


proximately  one  x-ray  diffraction  fringe.  Objects  1,  2,  and  3 
are  clearly  visible  and  are  marked  with  identical  shading, 
indicating  motion  in  the  same  direction  under  the  influence 
of  the  acoustic  radiation  pressure.  The  objects  identified  by 
arrows  in  Figs.  4(a)-4(c)  are  no  longer  visible  in  the  image 
and  can  therefore  be  identified  as  either  artifacts  arising  from 
optical  imperfections  in  the  sample  container  or  imaging  op¬ 
tics,  or,  as  objects  lying  outside  the  ultrasonic  field. 

The  salient  feature  of  the  acoustic  contrast  x-ray  imaging 
described  here  is  its  ability  to  highlight  density  variations: 
both  the  x-ray  phase  contrast  and  the  movement  of  the  se¬ 
lected  object  respond  to  density  variations.  Inherent  in  the 
method  is  the  rejection  of  all  features,  whether  absorptive  or 
refractive,  from  any  object  that  does  not  respond  to  ultra¬ 
sonic  radiation  pressure,  or  any  object  within  a  body  that  is 
not  within  the  ultrasonic  field.  This  feature  gives  the  diag¬ 
nostician  the  ability  to  address  a  specific  region  within  a 
body  with  directed  ultrasound  and  to  obtain  an  x-ray  image 
of  that  region  only.  Without  a  significant  reduction  in  the 
resolution  of  the  present  imaging  system  of  2.5  /xm/pixel  at 
the  object  plane,  a  field  of  view  of  10-15  cm  in  diameter  is 
possible  with  current  technology  by  using  a  larger  imaging 
screen.  The  present  results  suggest  application  of  the  method 
to  detection  of  small-sized  tumors  and  other  lesions  distin¬ 
guished  from  surrounding  tissue  by  their  elastic  properties 
and  density  differences  as  small  as  a  few  percent. 
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Abstract — The  current  in  a  circuit  produced  by  a  time- varying  polarization  generated  by  an  object  within  a  body 
in  response  to  the  passage  of  ultrasound  (US)  is  found  through  solution  to  Maxwell’s  equations.  Current 
generation  is  modelled  in  a  body  with  uniform  dielectric  constant  and  conductivity  where  a  finite  region  within 
the  body  produces  an  ultrasonic  vibration  potential.  The  results  of  the  calculation  form  the  mathematical  basis 
for  imaging  through  use  of  the  ultrasonic  vibration  potential.  (E-mail:  Gerald_Diebold@brown.edu)  ©  2005 
World  Federation  for  Ultrasound  in  Medicine  &  Biology. 
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INTRODUCTION 

The  generation  of  a  voltage  in  a  colloidal  or  ionic 
suspension  by  ultrasound  (US)  is  known  as  the  ultra¬ 
sonic  vibration  potential.  In  the  case  of  a  colloid, 
particles  suspended  within  solution  are  charged  and, 
in  turn,  are  surrounded  by  a  counter  charge  in  the 
fluid  giving  the  solution  overall  charge  neutrality. 
Because  the  density  of  the  particle  and  the  surround¬ 
ing  fluid  typically  differ  from  each  other,  the  pas¬ 
sage  of  a  sound  through  a  suspension  causes  the 
motions  of  the  particles  and  fluid  to  differ:  in  the 
typical  case  where  the  particles  are  denser  than  the 
surrounding  fluid,  the  particles  will  execute  a  smaller 
amplitude  motion  than  that  of  the  fluid,  with  a  phase 
that  lags  that  of  the  fluid.  The  flow  of  the  fluid  relative 
to  the  particle  distorts  the  normally  spherical  charge 
distribution  of  the  counter  charge  and  results  in  the 
creation  of  a  dipole  at  the  site  of  each  particle,  which 
alternates  in  direction  on  alternating  cycles  of  the 
sound  wave.  If  a  pair  of  electrodes  is  placed  in  solu¬ 
tion  a  distance  one  half  of  the  sound  wavelength  apart, 
an  alternating  voltage  can  be  recorded  across  the  elec¬ 
trodes.  The  theory  of  the  vibration  potential,  which 
was  originally  predicted  by  Debye  (1933)  for  ionic 
solutions,  has  been  reviewed  by  several  authors 
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(Hunter  1998;  Dukhin  and  Goetz  2001;  Zana  and 
Yeager  1983). 

Here,  we  discuss  the  generation  of  the  vibration 
potential  as  a  means  of  imaging  of  colloidal  or  ionic 
species  within  an  inert  medium  (Beveridge  et  al. 
2004a,  2004b).  The  problem  is  to  calculate  the  current 
developed  in  a  pair  of  electrodes  when  US  passes 
through  a  colloidal  or  ionic  object  within  a  body. 
The  calculation  for  the  Dirichelet  problem  is  carried 
out  for  a  pair  of  parallel  electrodes  at  zero  poten¬ 
tial  in  contact  with  a  medium  with  an  isotropic  con¬ 
ductivity  and  dielectric  constant,  within  which  is 
embedded  an  object  with  an  arbitrary  geometry  of 
the  same  conductivity  and  dielectric  constant  as  the 
surrounding  medium,  but  which  produces  a  vibration 
potential  on  the  passage  of  US.  There  are  two 
mechanisms  for  production  of  an  alternating  current 
in  the  circuit  connecting  the  two  electrodes:  the  first 
is  the  production  of  a  current  driven  by  the  electric 
field  in  the  body  induced  by  the  polarization,  and 
the  second  is  the  generation  of  charges  on  the  elec¬ 
trodes  as  a  result  of  imposition  of  a  time-varying 
polarization  in  the  region  between  the  parallel  plates. 
This  paper  is  divided  into  sections  giving  a  derivation 
of  Poisson’s  equation  for  determining  the  field  from  an 
arbitrary  colloidal  object,  derivation  of  a  Green’s 
function  for  the  parallel  electrode  geometry,  and  der¬ 
ivation  of  three  expressions  that  can  be  used  for  cal¬ 
culating  the  current  in  an  external  circuit  for  an  arbi¬ 
trary  object. 
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POISSON’S  EQUATION 

The  electric  field  E  produced  across  a  region  of 
colloid  by  a  pressure  p ,  according  to  O’Brien  (1988)  is 
given  by: 


E  = 


f^PmP-E 

PmCT* 


Vp, 


(1) 


where /is  the  volume  fraction  of  colloidal  particles  in  the 
fluid,  p,E  is  the  electrophoretic  mobility  of  the  particles, 
pm  is  the  mass  density  of  the  suspension  and  Apm  is  the 
difference  in  densities  between  the  particles  and  the 
surrounding  fluid,  and  <x*  is  the  complex  conductivity 
(O’Brien  1986;  The  complex  conductivity  <x*  used  by 
O’Brien  is  a  generalized  conductivity  that  relates  the 
current  density  to  the  electric  field  in  a  colloid). 

Consider  generation  of  the  vibration  potential  in  a 
medium  with  conductivity  cr,  and  dielectric  constant  e, 
assumed  uniform  throughout  the  medium  under  consid¬ 
eration,  but  with  one  region  of  the  medium,  referred  to 
here  as  the  object,  that  generates  a  vibration  potential  on 
the  passage  of  sound.  By  combining  the  charge  conser¬ 
vation  relation  v  •  J  =  dp/dt ,  where  J  is  the  current 
density  and  p  is  the  charge  density,  with  the  constitutive 
relation  J  =  crE,  where  E  is  the  electric  field  and  a  is  the 
real  part  of  the  conductivity,  and  Maxwell’s  equation  p 
=  V  •  D,  where  D  is  the  electrical  displacement,  it  fol¬ 
lows  that: 


V-(oE)=-^(V-D).  (2) 


V- 


oE  +  —  (sE  +  Pa) 


=  0. 


(5) 


For  the  present  problem,  where  the  media  under 
consideration  would  normally  be  biologic  tissue,  the 
frequencies  where  the  vibration  potential  can  be  gener¬ 
ated  are  restricted  owing  to  ultrasonic  absorption;  hence, 
it  is  reasonable  to  make  the  quasi-electrostatic  approxi¬ 
mation  where  V  X  E  —  0,  so  that  the  electric  field  can  be 
written  in  terms  of  a  potential  4>  as: 


E=  —  V(f>. 


(6) 


Substitution  of  eqn  (6)  into  eqn  (5)  gives: 


^+4,v2*=r,<vp-)- 


m 


from  which  the  potential  in  the  time  domain  can  be 
determined  from  a  knowledge  of  Pa  in  space  and  time. 
Equation  (7)  can  be  found  in  the  frequency  domain  by 
assuming  that  the  acoustic  pressure,  the  polarization  and 
the  potential  vary  sinusoidally  in  time  according  to: 

p  =  Re[p(x)ei(0t] 


Pfl  =  Re[Pa(x>>1 


c j>  =  Re[(Kxyw], 


which,  on  substitution  into  eqn  (7),  gives  Poisson’s  equa¬ 
tion  for  the  frequency  domain  potential  in  terms  of  the 
acoustic  polarization  as: 


The  effect  of  the  passage  of  sound  through  the 
object,  which  is  taken  to  be  colloidal,  is  to  produce  an 
acoustic  polarization  Pa  within  the  object,  so  that  the 
displacement  vector  becomes: 

D  =  sE  +  Pa,  (3) 

where  s  is  the  dielectric  constant  of  the  object  and  the 
surrounding  medium,  assumed  to  be  identical.  It  is  as¬ 
sumed  that  the  polarization  P  varies  in  space,  depending 
on  whether  the  coordinate  is  within  the  object  or  in  the 
inert  region  of  the  body.  The  polarization  produced  by 
the  sound,  within  the  context  of  the  present  calculation, 
is  taken  to  be: 


P a  =  a(x,M)Wp  (4) 


where 


V24>  = 


id) 

cr  +  icos 


V  P 


a’ 


(8) 


or 


V24>=  V  •(atVp),  (9) 

where  ot  is  given  by: 

.  /ft) 

a/x,^))  =  — — — a(x,  co).  (10) 

cr  +  nos 

It  can  be  seen  that  eqns  (8)  and  (9)  are  Poisson’s  equa¬ 
tions  with  the  acoustic  polarization  acting  as  a  source. 
The  source  is  dependent  on  derivatives  of  both  the  col¬ 
loidal  object  and  the  pressure  in  space. 


GREEN’S  FUNCTION  FOR  THE  POTENTIAL 


«(x,  w)  =  sf  \pmp,E/ pmcr*  . 

The  notation  for  ai(x,ft>)  explicitly  contains  the  space 
variable  to  indicate  that  it  takes  on  different  values  in  the 
colloidal  and  non  colloidal  regions.  The  combination  of 
eqn  (3)  with  eqn  (2)  gives: 


Consider  the  geometry  for  imaging  a  colloidal  ob¬ 
ject  of  arbitrary  shape  inside  a  weakly  conducting  body 
with  a  pair  of  conducting  electrodes  in  contact  with  the 
external  surface  of  the  body,  as  shown  in  Fig.  1.  The 
electrodes  are  considered  to  be  parallel  and  spaced  a 
distance  h  apart,  with  the  direction  of  the  US  taken  to  be 
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Fig.  1.  A  colloidal  object  is  embedded  within  a  body.  The 
object  and  the  surrounding  region  are  taken  to  have  a  conduc¬ 
tivity  a  and  a  dielectric  constant  e.  Only  the  colloidal  region 
within  the  body  produces  a  vibration  potential;  the  region 
outside  the  object  is  considered  to  be  inert. 


along  the  z-axis.  The  electrodes  are  assumed  to  be  con¬ 
nected  to  each  other  through  a  conductor  or,  equiva¬ 
lently,  through  an  amplifier  with  an  input  impedance  of 
nearly  zero  so  that  the  potentials  of  each  of  the  electrodes 
are  effectively  zero.  The  potential  c/>  within  the  body  can 
be  determined  from  a  Green’s  function  solution  to  Pois¬ 
son’s  equation  through  use  of  a  Dirichlet  Green’s  func¬ 
tion  Gd(x,x'),  that  satisfies  the  equation: 

V2Gd(x,  x')  =  8(x  -  x') ,  (11) 

where  x  and  x'  are  the  field  and  source  points,  respec¬ 
tively.  The  potential,  according  to  Green’s  theorem  is 
then  given  by: 


<?i(x)  =  J  GdV  ■  (a'^'p')dY  +  § 


4> 


YGr 


dri 


t~Gd 


d4> 

dn' 


dS\ 

(12) 


where  the  integration  in  both  integrals  is  over  the  primed 
coordinates,  and  where  the  derivatives  in  the  surface 
integral  are  taken  with  respect  to  the  surface  normal  n. 
Because  the  potentials  on  the  planes  at  z  =  0  and  z  =  h 
are  taken  to  be  zero,  it  is  convenient  to  construct  a 
Green’s  function  that  is  also  zero  on  the  surfaces  of  the 
electrodes  and  to  take  the  surface  of  integration  to  be  a 
rectangular  enclosure  in  the  body,  with  its  ends  at  the 
positions  of  the  electrodes  where  the  potential  is  zero, 
and  where  the  enclosure  extends  to  infinity  in  the  x  and 
y  directions. 

Consider  the  Fourier  transformation  of  the  spatial 
coordinate  in  eqn  (11)  through  the  transform  pair: 


f(x,y)  =  ID-  (13b) 

—  00  —00 


to  give  the  Green’s  function  Gd(kx ,  Ky ,  z\  x')  in  reciprocal 
space.  Fourier  transformation  of  eqn  (11)  gives: 


Ky,  z;  x') 


1 

(277? 


S(z  ~  &e-i(**'+K’y0, 


(14) 


where  the  Green’s  function  in  reciprocal  xy  space  is 
Gd(kx,  Ky,  z).  Because  the  Green’s  function  must  vanish 
at  z  =  0  and  z  =  h,  it  can  be  of  the  form  GD  = 
2Amsin(dmz),  where  dm  =  mir/h,  Am  is  a  constant  inde¬ 
pendent  of  z  and  m  is  an  integer  ranging  from  1  to  o°.  By 
multiplication  of  eqn  (14)  by  sin  dm>z  and  integration  of 
the  resulting  expression  from  0  to  h,  it  can  be  shown  that 
the  reciprocal  space  Green’s  function  is  given  by: 


Gd(kx,  Ky,  z;  x') : 


27T2h 


/ rmr  \  / mu 

sinlxz)sinllTz/ 


/  rmr 

(x 


I  -  1  +  K?  +  kI 


x  e-iK*x'-iKyy\  (15) 


which,  on  transformation  back  to  coordinate  space  ac¬ 
cording  to  eqns  (13a)  and  (13b)  gives  the  Green’s  func¬ 
tion  as: 


Gl)(x,  x')  =  -  27J2h  2  sin 

1 


mir  \ 

h  zr 


x 


giKx(x-x')+iKy(y-y') 


+  Kr  +  Kv 


dKxdKy. 


(16) 


Because  the  potential  is  zero  on  the  upper  and 
lower  electrodes  by  assumption,  and  Gd(x,  x')  is  zero  on 
the  upper  and  lower  electrodes  by  construction,  the  sur¬ 
face  integrals  in  eqn  (12)  are  zero,  so  that  the  frequency 
domain  potential  is  given  by  the  volume  integral: 


4>(X)  =  Jr  Gd(x,  x')V  ■  [a'(x',  0j)Vp(x’)]dV,  (17) 


00  00 

f(Kx,  Ky)  =  (2?J  /  f(x,  y)e-'(K^dxdy  (13a) 


where  Gd(x,x’)  is  given  by  eqn  (16)  and  where  the 
integration  over  the  primed  coordinates  is  explicitly  in¬ 
dicated  in  the  integrand. 
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THE  CURRENT  IN  THE  ELECTRODES 

There  are  two  mechanisms  for  production  of  current 
in  the  electrodes  (O’Brien  1988).  If  the  region  between 
the  two  plates  is  a  nonconducting  dielectric,  a  current 
will  be  generated  in  the  external  circuit  connecting  the 
two  electrodes  as  a  result  of  the  flow  of  charge  induced 
by  the  time- varying  polarization  induced  by  the  US  in 
the  colloidal  or  ionic  object.  The  current  density  JP 
generated  by  a  time-varying  displacement  vector  D  is 
given  by  JP  =  dD/dt,  where  the  displacement  vector  can 
be  found  from  the  potential  through  eqn  6  and  the  con¬ 
stitutive  relation  for  a  dielectric  D  =  sE.  Thus,  JP  is 
found  from  the  potential  as 


J p(x)=-eytV<Kx).  (18) 

Similarly,  when  a  field  is  produced  in  a  weak  con¬ 
ductor  by  a  time- varying  polarization,  the  current  density 
for  the  free  charge  JF  is  related  to  the  electric  field  as  JF 
=  crE,  so  that  the  current  density  expressed  in  terms  of 
the  potential  becomes: 

Jf(x)  =  —  <tV  (j)(x),  (19) 

The  current  produced  by  both  mechanisms  I  is  found  by 
integrating  the  normal  component  of  the  electric  field, 
that  is,  the  component  perpendicular  to  one  of  surfaces, 
over  the  surface  of  one  of  the  electrodes  5;  hence,  the 
expression  for  the  frequency  domain  current  /  in  the 
external  circuit  is  given  by: 


Note  that  the  factor  of  (cr  +  ia>s)  appears  both  in  the 
prefactor  to  the  integral  and  in  the  denominator  of  eqn 
(10)  that  defines  o';  hence,  eqn  (20)  can  be  expressed  in 
terms  of  a  rather  than  o' .  The  integration  over  x  and  y 
gives  (2tt)2S(kx)8(k)9  which  simplifies  eqn  (20)  consid¬ 
erably.  Through  use  of  the  identity: 


sm nx 

2— (-i)"  = 


eqn  20a  can  be  shown  to  reduce  to: 


/=-tz 


J‘ 


dxfdf 


d  /  d  _ 

+  xx[  a—p 
dy\  dy 


d  /  d  _ 
+  dAadzP 


dxdydz ,  (21) 


where  the  primes  have  been  dropped  in  the  integrand.  If 
each  of  the  terms  in  eqn  (21)  is  integrated  by  parts, 
noting  that  the  dimensions  of  the  object  are  such  that  it 
does  not  contact  either  of  the  electrodes,  so  that  a  can  be 
taken  as  zero  at  z  =  h  or  z  =  0,  it  is  not  difficult  to  show 
that  the  only  term  that  remains  is  the  third  term  in  the 
integrand,  so  that  the  current  is  given  by: 


ico  p  8 

I  =  -jj- J  v  a(x ,  co)—  p(x)dxdydz.  (22) 


It  is  also  possible  to  carry  out  an  integration  of  this 
expression  by  parts  as  well,  giving  an  alternate  expres¬ 
sion  for  the  current  in  the  circuit  as: 


/  = 


(cr  +  i(oe) 


f  f  dxdyfV'FGD(x,x')\z=0 


io)  r  _ 

I  =  — ^  J  v  [V a(x,  co)]^  p(x)dxdydz.  (23) 


V' 

•  [at(x/,  a))Yp(x')]d 

V\ 

which,  on  substitution  of  Gd(x,  x  )  from  eqn  (16)  carry¬ 
ing  out  the  differentiation  with  respect  to  z,  and  evalua¬ 
tion  of  the  integral  at  z  =  0  gives: 


A  useful  integral  can  be  derived  from  eqn  (22)  by 
using  the  linearized  Navier-Stokes  law  for  an  inviscid 
fluid,  Vp/p  =  —  u  where  u  is  the  displacement  of  the 
medium  in  the  acoustic  field.  The  acceleration  in  a  plane 
wave  can  be  expressed  in  terms  of  the  displacement 
through  the  relation  ii  =  —  &>2u.  Thus,  an  expression  for 
the  current  in  terms  of  the  displacement  in  the  plane 
wave  can  be  written  as: 


(cr  +  io)s)  r  r  r  JL  /  rmr\  ( mir 


X 


JJr 


JKx(x-x')+iKy(y-y') 


miry 

7T 


^dKxdKyV 


+  +  Ky 


•  [at(x/,  co)V'/?(x')]c/ 

1/.  (20) 


ip  of  p  _ 

/  =  — —  J  v  a(x ,  of)u(x)dxdydz.  (24) 

INFINITE  SLAB 

The  problem  of  a  slab  of  thickness  a  located  be¬ 
tween  two  metal  electrodes  connected  by  a  low-imped¬ 
ance  amplifier,  as  shown  in  Fig.  2,  can  be  solved  by 
application  of  eqn  (23).  The  impedance  of  the  amplifier 
is  taken  to  be  so  small  that  it  acts  as  a  short  circuit,  so  that 
the  potential  of  the  two  electrodes  can  be  considered  as 
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electrode  cp=0 

t>~ 

preamplifier 

c 

a 

colloidal  layer 

b 

electrode  cp=0 

y 


m  -  m  =  0  =  - f  Edz  =  —\  f  Dgdz 

J0  s  Jo 

/>a+b 

P  adz 

h 


(27) 


Therefore,  the  component  of  the  displacement  vector 
along  the  z-axis  in  all  three  layers  is  given  by: 


1  rb 

D°=hSb 


P adz. 


(28) 


Fig.  2.  An  infinite  slab  of  a  colloidal  suspension  lies  between 
two  electrodes.  The  input  impedance  of  the  amplifier  is  con¬ 
sidered  to  be  sufficiently  low  that  the  potential  of  each  elec¬ 
trode  can  be  taken  to  be  zero. 


zero.  The  slab  and  the  surrounding  material  are  consid¬ 
ered  to  have  the  same  dielectric  constant  s;  however,  the 
slab  in  region  “a”,  located  between  z  =  b  and  z  =  b  + 
a  is  assumed  to  be  the  only  region  where  a  vibration 
potential  is  generated.  The  ultrasonic  beam  is  assumed  to 
be  continuous  and  to  propagate  along  the  positive  z-axis. 
For  the  present  problem,  the  distribution  of  colloid  in 
space  is  given  by  o:(z,m)  =  a0[u(z  —  b)  -  u  (z  —  b  —  a)], 
where  a0  refers  to  a  value  of  o:(x,m)  independent  of  x  in 
layer  a  whose  magnitude  is  given  by  eqn  (4),  and  u  is  the 
Heaviside  function.  Because  the  layers  are  considered  to 
be  infinite  in  the  x  and  y  directions,  a  current  per  unit  area 
J  is  generated  by  the  time- varying  polarization,  which, 
according  to  eqn  (23),  is  found  to  be: 

i(OGLo 

J= - +  (25) 

It  is  possible  to  consider  the  same  problem  without 
the  use  of  eqn  (23)  by  assuming  a  time- varying  polar¬ 
ization  Pa  that  is  generated  by  the  vibration  potential,  and 
calculating  the  resulting  current  when  the  material  in 
layers  “b”  and  “c”  are  dielectrics  with  dielectric  constant 
s.  From  eqn  (3  a),  the  electric  field  in  the  three  regions  is 
given  by: 

E(z)  =  '[A,  -  Pa(z)],  (26) 

where  E ,  D  and  Pa  are  the  components  of  electric  field, 
displacement,  and  acoustic  polarization  along  the  posi¬ 
tive  z-axis.  The  acoustic  polarization  has  a  value  Pa 
inside  region  a ,  and  is  zero  in  regions  b  and  c.  Because 
there  are  no  free  charges  in  the  three  layers,  the  displace¬ 
ment  vector  has  a  constant  value  D0.  From  eqn  (6),  it 
follows  that: 


It  is  not  difficult  to  show  that  the  electric  fields  and 
polarizations  in  regions  b  and  c  are  identical.  If  an 
integration  volume  is  taken  as  a  closed  cylindrical  vol¬ 
ume  crossing  the  region  between  the  upper  electrode  and 
region  b ,  then  it  follows  from  Maxwell’s  equation  v  •  D 
=  p  that  D0=  —  c ts ,  where  as  is  the  charge  density  on  the 
upper  electrode.  Because  magnitude  of  the  current  den¬ 
sity  J  is  given  by  /  =  dajdt ,  the  frequency  domain 
current  density  is: 

lOJ  f*b+a 

J=~h\b  P^,  (29) 

which,  for  an  acoustic  polarization  Pa  =  a(co,  x)Vp, 
recovers  eqn  (25),  as  found  directly  from  eqn  (23), 
above. 

Note  the  simplification  that  follows  from  the  as¬ 
sumption  of  a  1-D  geometry  in  this  problem.  If  an 
acoustic  beam  with  a  finite  cross-section  were  employed, 
then  the  direction  of  the  electric  field  would  not  be 
parallel  to  the  z-axis  over  the  entire  area  of  the  slab,  and 
a  more  complicated  calculation  of  the  current  density 
would  be  required.  The  magnitude  of  the  error  in  deter¬ 
mining  the  current  density  through  use  of  eqn  (25)  for  a 
finite  cross-section  acoustic  beam  would,  of  course,  de¬ 
pend  on  the  dimensions  of  the  beam  relative  to  the  width 
of  the  gaps  between  the  colloidal  layer  and  the  conduc¬ 
tors. 

DISCUSSION 

Equations  (22)  to  (24)  give  three  equivalent  expres¬ 
sions  for  the  current  in  a  circuit  generated  by  a  colloidal 
or  ionic  object  inside  a  body  irradiated  by  a  beam  of  US. 
The  third  of  these  gives  perhaps  the  most  intuitive  pic¬ 
ture  of  the  process  of  current  generation,  in  that  it  states 
that  the  current  is  proportional  to  material  displacement 
in  the  acoustic  wave  integrated  over  the  region  where 
colloid  is  present.  Thus,  the  production  of  current  is  seen 
as  proportional  to  the  integration  of  the  dipole  moments 
at  the  sites  of  the  particles  over  the  region  where  there  is 
a  displacement  caused  by  the  acoustic  wave.  Note  that, 
depending  on  the  frequency  and  spatial  extent  of  the 
wave,  the  integration  can  be  finite  or  zero.  This  property 
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is  seen  explicitly  in  the  solution  for  the  current  given  by 
eqn  (25),  which  shows  that  a  sinusoidal  pressure  wave 
whose  spatial  extent  is  smaller  than  the  length  of  the 
colloidal  region  generates  no  signal  when  the  wave  is 
confined  within  the  colloidal  layer  (in  which  case  p(b  + 
a)  =  p(b)  =  0);  only  when  such  a  short  burst  enters  or 
exits  the  layer  is  a  current  produced.  For  a  continuous 
wave,  or  an  ultrasonic  burst  whose  spatial  extent  is  large 
compared  with  the  layer  thickness,  when  an  even  or  odd 
number  of  half  wavelengths  of  the  acoustic  wave  equals 
the  length  of  the  colloidal  layer,  the  vibration  potential 
across  the  layer  integrates  to  zero,  or  takes  on  a  maximal 
value.  For  other  wavelengths,  the  signal  integrates  to 
some  current  intermediate  between  these  two  extrema. 
Further,  it  is  easy  to  see  that  recording  the  frequencies 
that  give  consecutive  maxima  and  minima  in  the  vibra¬ 
tion  potential  permits  determination  of  the  layer  thick¬ 
ness.  The  simple  example  of  the  slab  also  makes  clear  the 
interpretation  of  eqn  (23),  in  that  a  gradient  of  the  con¬ 
centration  of  colloid  is  required  for  a  current  to  be 
generated  for  burst  of  US:  a  fixed  amplitude  sound  wave 
crossing  a  boundary  between  colloidal  and  inert  material 
generates  a  current  at  the  points  where  it  enters  or  leaves 
the  colloidal  region.  Equally,  eqn  (22)  states  that,  if  there 
is  an  acoustic  wave  that  varies  in  space  propagating 
wholly  within  a  colloidal  layer  of  uniform  concentration, 
it  is  possible  to  observe  a  current  from  the  vibration 
potential. 

Different  methods  of  generating  an  image  are  evi¬ 
dent  in  the  results  obtained  here.  In  addition  to  determi¬ 
nation  of  the  thickness  of  a  colloidal  or  ionic  layer  in  a 
1-D  geometry  already  mentioned,  a  method  based  on  the 
time  of  appearance  of  the  vibration  potential  is  inherent 
in  eqns  (23)  and  (24).  The  distance  a  colloidal  region  lies 
from  the  launching  point  of  an  ultrasonic  wave  can  be 
determined  from  a  knowledge  of  the  sound  speed  in  the 
body  and  the  appearance  time  of  the  vibration  potential 
signal  relative  to  the  launching  time  of  the  sound  wave, 
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in  a  manner  analogous  to  conventional  ultrasonic  imag¬ 
ing.  In  general,  the  frequency-dependence  of  the  current 
for  a  continuous  ultrasonic  wave,  or  the  time-dependence 
of  any  acoustic  burst  sequence  generated  by  an  object 
whose  geometry  is  known,  can  be  determined  from  eqn 
(23).  It  is  expected  that,  for  a  number  of  simple  geome¬ 
tries,  straightforward  expressions  for  the  vibration  poten¬ 
tial  can  be  obtained  in  closed  from,  which  will  prove  of 
value  in  the  interpretation  of  vibration  potential  signals 
for  imaging  (Beveridge  et  al.  2004b).  In  summary,  the 
present  work  gives  explicit  expressions  for  the  current 
recorded  in  an  imaging  experiment  for  objects  whose 
geometry  is  known  a  priori ;  in  addition,  it  gives  a  start¬ 
ing  point  for  solution  of  the  inverse  problem  where 
vibration  potential  data  are  gathered  to  determine  the 
spatial  extent  of  an  arbitrary  object  whose  geometry  and 
location  within  a  body  is  sought. 
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Abstract 

We  report  the  use  of  ultrasonic  radiation  pressure  with  phase  contrast  x-ray 
imaging  to  give  an  image  proportional  to  the  space  derivative  of  a  conventional 
phase  contrast  image  in  the  direction  of  propagation  of  an  ultrasonic  beam. 
Intense  ultrasound  is  used  to  exert  forces  on  objects  within  a  body  giving 
displacements  of  the  order  of  tens  to  hundreds  of  microns.  Subtraction  of 
images  made  with  and  without  the  ultrasound  field  gives  an  image  that  removes 
low  spatial  frequency  features  and  highlights  high  frequency  features.  The 
method  acts  as  an  acoustic  ‘contrast  agent’  for  phase  contrast  x-ray  imaging, 
which  in  soft  tissue  acts  to  highlight  small  density  changes. 


1.  Introduction 

When  an  x-ray  beam  traverses  a  body  the  spatial  intensity  profile  of  the  beam  is  modified  by 
the  effects  of  absorption,  scattering  and  phase  changes,  all  of  which  depend  essentially  on  the 
density  of  electrons.  The  first  two  of  these  effects  form  the  basis  of  conventional  radiography 
where  a  shadow  is  produced  in  the  image  wherever  an  object  is  of  higher  density  than  its 
surroundings.  A  second  contrast  mechanism,  known  as  phase  contrast,  follows  directly  from 
the  Huygens-Fresnel  plane-wave  theory  of  image  formation,  which  responds  to  changes  in 
index  of  refraction.  The  effect  of  phase  contrast  is  subtle  and  is  seen  only  when  a  spatially 
coherent  source  of  x-rays  is  employed.  The  Fresnel-Kirchhoff  integral  (Born  and  Wolf  1980, 
Cowley  1984)  for  one  dimension  in  the  image  plane  x,  gives  the  field  amplitude  f(x)  of  the 
x-ray  beam  at  a  distance  z  from  the  object  as 

exp(— i kz)  J  q(X)  exp  ^  ~  j  dX  (1) 

where  y  is  the  wavelength  of  the  radiation,  k  is  the  wave  number  of  the  radiation  and  q  ( X ) 
describes  the  transmission  of  the  object.  The  transmission  function,  q(x),  for  an  object 
gives  the  effect  of  the  object  on  both  the  phase  and  the  amplitude  of  the  incoming  wave 
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and  determines  the  intensity  pattern  in  the  image.  For  a  thin  object  with  a  small  integrated 
absorption,  q{x)  can  be  written  as 

q(x)  =  exp(i </>(x)  -  /x(*))  (2) 

where  <p  and  /x  describe  the  change  in  phase  and  amplitude  of  the  incoming  wave,  respectively, 
in  passing  through  the  object.  The  intensity  of  the  image  I(x)  =  \f{x)\2  for  a  weak  phase 
object  can  be  shown  to  be 

/(*)  =  !-  (3) 

In 

Similarly,  for  a  weak  absorber,  the  intensity  at  the  image  plane  is  given  by 

I(x)  =  1  —  2  /x(v).  (4) 

In  a  typical  phase  contrast  image  of  an  object  embedded  in  a  lower  density  medium,  the 
usual  shadowgraph  from  absorption  is  present  according  to  equation  (4),  but  at  the  edges  of 
the  object  interference  fringes  are  observed.  According  to  equation  (3),  which  is  valid  for 
low  spatial  frequencies,  the  intensity  in  the  image  varies  as  the  second  derivative  of  the  phase 
variation  in  the  object,  which  tends  to  be  greatest  for  smooth  objects  at  the  interface  between 
the  object  and  the  surrounding  medium. 

Several  techniques  have  been  used  to  record  phase  variations.  In  electron  microscopy, 
for  example,  phase  contrast  is  obtained  by  defocusing  the  electron  beam  and  is  imaged 
on  a  scintillation  screen.  X-ray  sources  with  high  spatial  coherence,  such  as  synchrotrons 
(Snigirev  et  al  1995,  Fulvia  Arfelli  et  al  2000,  Momose  1995),  microfocus  x-ray  tubes 
(Pogany  et  al  1997,  Wilkins  et  al  1996)  or  laser  plasma  x-ray  sources  (Krol  et  al  1997a,  1997b) 
are  commonly  used  for  phase  contrast  imaging. 

A  detailed  analysis  of  the  relative  contrast  produced  in  an  image  for  a  fixed  difference 
in  density  shows  that  the  phase  contrast  is  more  sensitive  than  the  absorption  for  detection  of 
density  differences  throughout  most  of  the  spectral  region  commonly  used  for  diagnostic  x-ray 
imaging  (Beckmann  et  al  1997).  For  instance,  the  thickness  of  water  needed  to  produce  a  1% 
phase  contrast  at  36  keV  photon  energy  is  2500  times  smaller  than  that  required  to  produce 
the  same  absorption  contrast  (Beckmann  et  al  1997).  As  a  consequence,  x-radiation  with  a 
high  photon  energy,  which  is  only  weakly  absorbed,  can  be  used  in  phase  contrast  imaging 
to  obtain  the  same  sensitivity  to  density  change  as  much  more  strongly  absorbed,  low  energy 
x-radiation  used  in  absorption  imaging,  the  obvious  advantage  of  the  former  being  a  far  lower 
dose  of  absorbed  radiation  in  the  irradiated  material.  The  size  of  the  smallest  phase  object  (i.e. 
an  object  where  /x  =  0)  that  can  be  detected  by  phase  contrast  can  be  estimated  for  an  object, 
with  index  of  refraction  ^object,  embedded  in  a  surrounding  medium,  with  index  of  refraction 
^medium,  from  the  relation 

4^1%  —  "  (^object  ^medium)  (5) 

A 

where  t\%  is  the  thickness  of  the  object  required  to  produce  an  intensity  change  of  1%, 
corresponding  to  a  phase  change  of  0.14  rad.  Consider  application  of  equation  (5)  to  tumour 
detection  in  soft  tissue.  Several  studies  (Johns  and  Yaffe  1987,  Duck  1990)  have  determined 
that  mass  density  differences  between  malignant,  for  instance  infiltrating  carcinomas,  and 
healthy  breast  tissue  are  of  the  order  of  10%  (Duck  1990,  Johns  and  Yaffe  1987).  Since  8  is 
approximately  10-6  for  soft  tissue,  the  quantity  <50bject  —  ^medium  is  of  the  order  of  10-7.  At 
an  x-ray  energy  of  40  keV,  corresponding  to  a  wavelength  of  3 1  pm,  equation  (5)  gives  t\%  as 
7  /xm,  which  implies  that,  ignoring  specific  morphological  features,  malignant  tissue  with  a 
thickness  of  about  10  /xm  can  be  distinguished  from  benign  tissue  assuming  a  signal-to-noise 
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ratio  of  at  least  100.  For  the  results  discussed  in  this  paper,  the  signal-to-noise  ratio  is  15, 
which  is  sufficient  to  distinguish  object  with  a  thickness  of  at  least  45  pan. 

Here  we  report  experiments  where  pulsed  ultrasound  is  used  to  create  displacements  of 
objects  within  water  and  soft  tissue.  The  displacement  of  objects  with  densities  different  from 
that  of  the  surrounding  body  is  used  to  modify  a  phase  contrast  x-ray  giving  an  image  with 
an  enhanced  phase  contrast  component.  Section  2  gives  calculations  of  the  radiation  pressure 
on  spheres  useful  for  predicting  displacements  of  spherical  objects  by  ultrasound.  Section  3 
describes  the  x-ray  apparatus  and  experiments  where  displacements  of  the  order  of  microns 
were  produced  by  acoustic  radiation  pressure  (Westervelt  1951,  Chu  and  Apfel  1982,  Beyer 
1978)  and  measured  phase  contrast  radiographically.  Section  4  describes  the  x-ray  acoustic 
method  of  image  formation  and  gives  examples  of  its  use.  The  discussion  section  comments 
on  some  of  the  features  of  the  method  and  discusses  applications. 


2.  Acoustic  radiation  pressure 


Acoustic  radiation  forces  arise  from  two  processes,  either  reflection  of  sound  by  an  object  as 
a  result  of  an  acoustic  impedance  change  (i.e.,  a  variation  in  either  density  or  sound  speed)  or 
by  absorption  of  sound  by  the  object.  A  number  of  studies  (Muthupillai  et  al  1995,  Gao  et  al 
1996,  Fatemi  and  Greenleaf  1998,  Nightingale  et  al  1999,  Sarvazyan  1998,  McAleavey  2003) 
have  shown  that  the  soft  tissue  can  be  manipulated  by  ultrasound  pressure.  Typically,  a 
mechanical  actuator  or  a  beam  of  ultrasound  moves  an  object  whose  position  is  monitored  in 
time  with  conventional  pulse-echo  ultrasonic  imaging.  It  has  been  shown  that  variations  in 
Young’s  modulus,  which,  in  part,  determines  sound  speed,  permit  acoustic  differentiation  of 
tissue.  It  is  known,  for  instance,  that  breast  tumours  have  a  significantly  different  Young’s 
modulus  from  surrounding  tissue,  and  that  tumours  move  as  rigid  bodies  in  response  to  acoustic 
radiation  pressure  (Gao  et  al  1996,  Nightingale  et  al  1999,  2002,  Sarvazyan  1998,  McAleavey 
2003).  It  should  also  be  noted  that,  at  least  in  a  mammographic  application,  the  collective 
motion  associated  with  respiration  and  blood  pressure  is  minimized  by  immobilizing  the 
breast. 

For  the  purposes  of  the  present  study,  where  reflection  of  the  acoustic  wave  is  considered, 
it  can  be  said  that  for  a  fixed  acoustic  impedance  change,  the  force  exerted  by  an  ultrasonic 
wave  is  proportional  to  the  time  average  of  the  energy  density  of  the  wave  and  the  area  of 
the  object  presented  to  the  sound  field.  An  expression  for  the  acoustic  radiation  force  on  an 
arbitrary  object  has  been  given  by  Westervelt  (1951)  and  evaluated  for  spheres  small  compared 
with  the  wavelength  of  the  radiation.  For  spheres  of  any  diameter,  King  (1934)  has  given  an 
expression  that  agrees  with  the  limiting  expression  given  by  Westervelt.  According  to  King, 
the  average  acoustic  radiation  force  exerted  on  a  sphere  of  radius  a  is  given  by: 


F 


ac 


2jtp2 

( ka)2k2c2p 


A(ka) 


(6) 


where  p  and  c  are  the  density  and  the  speed  of  sound  of  the  medium,  respectively,  p  is  the 
acoustic  pressure,  k  is  the  wave  number  of  the  radiation,  and  A(ka )  is  given  by 


A(ka)  = 


1 


2 


[(ka)2~  3(1  ~£)]2 

(kaf 


+  E 


/  (n  +  1)  (( ka )2  —  n(n  +  1))2\ 

\HnH„+t  (£a)4(n+1)  ) 
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where  the  functions  Hn  for  small  ka  are  given  by 


n  —  0: 


Ho 


(1  +  (ka)2) 
( ka )2 


n  —  1: 


n  >  2: 


Hn  = 


(»  +  D2(n”=i  (2«  -  D)2 


(fctf) 


4«+2 


1  + 


77  —  1 


(n  +  \)(2n  -  1) 


(ka)2  +  •  • 


where  p  is  the  density  of  the  fluid  and  ps  is  the  density  of  the  sphere.  For  ka  >  2,  the  functions 
Hn  are  given  by 


n  =  0: 

(1  +  (ka)2) 
0  (ka)2 

n  —  1: 


where  Jn  is  a  Bessel  function  of  order  n,  and  where  the  argument  of  each  of  the  Bessel  functions 
is  ka. 

Figure  1  is  a  plot  of  force  versus  ka  from  equation  (7),  with  the  inset  showing  a  portion  of 
the  curve  for  small  ka.  The  parameters  used  in  equation  (7)  were  for  a  Teflon  bead  with  a  = 
1.5  mm,  p/ ps  =  2;  and  for  water  with  c  =  1500  m  s-1  and  p  =  103  kg  m-3.  The  curves  in 
figure  1  show,  unlike  optical  radiation  forces,  a  steady  increase  in  the  force  with  increasing 
ka;  the  acoustic  radiation  pressure,  of  course,  has  a  much  weaker  dependence  on  ka. 


3.  Radiation  force  experiments 

Experiments  were  carried  out  to  determine  the  deflection  of  a  Teflon  bead  in  water  using 
an  x-ray  imaging  apparatus  to  determine  the  deflection  of  the  bead  as  a  function  of  acoustic 
pressure.  The  x-rays,  generated  by  a  microfocus  x-ray  tube  (Oxford  Ultrabright  Microfocus 
UB-M1),  were  directed  onto  a  sample  cell  consisting  of  a  3  cm  diameter  PVC  tube  which  had 
portions  machined  out  and  replaced  with  Mylar  foil  to  reduce  x-ray  intensity  losses.  The  x-ray 
beam,  as  shown  in  figure  2,  propagated  through  a  He  filled  tube  to  a  Gd202S(Tb)  fibre  optic 
scintillation  plate  (Hamamatsu,  Inc.  Model  J6676)  the  fluorescence  from  which  was  imaged 
onto  a  liquid  nitrogen  cooled  CCD  camera  (Roper  Scientific  Model  7382-0001).  The  intensity 
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Figure  1.  Radiation  force  (in  N)  versus  ka  (dimensionless)  from  equation  (6).  Inset:  magnified 
view  of  the  first  part  of  the  curve.  The  portions  of  the  curves  for  ka  <2  were  calculated  using  the 
small  ka  expressions  for  the  functions  Hn. 


Figure  2.  Diagram  of  the  experimental  apparatus.  X-radiation  generated  by  a  microfocus  tube 
penetrates  a  sample  and  is  detected  by  a  CCD  camera  that  views  a  scintillation  plate.  The  CCD 
camera  is  read  by  the  computer  which  stores  the  images  and  performs  the  subtractions. 


maximum  of  the  emission  x-ray  spectrum  was  40  keV.  The  source- to -detector  distance  was 
2.6  m,  while  the  sample-to-source  distance  was  0.35  m,  giving  a  magnification  of  7.5  and  a 
maximum  phase  contrast  (Pogany  et  al  1997)  for  spatial  frequencies  corresponding  to  4.3  /xm. 

The  acoustic  transducer  used  to  produce  radiation  pressure  was  a  1  MHz,  LiNbC>3- 
transducer  located  3  cm  from  the  bead.  One  end  of  the  PVC  sample  tube  was  terminated 
by  a  synthetic,  sound-absorbing  plastic.  The  pulse  train  from  the  function  generator  (Agilent 
Model  33250A)  used  to  drive  the  transducer  consisted  of  15  /i s  bursts  at  a  repetition  rate  of 
133  Hz.  The  signal  from  the  function  generator  was  amplified  by  a  power  amplifier  with  a 
peak  output  power  of  1.5  kW  and  delivered  to  the  transducer.  Although  the  low  duty  cycle  of 
the  pulse  sequence,  approximately  0.2%,  resulted  in  a  reduced  time  average  acoustic  power 
from  the  transducer,  the  resulting  radiation  pressure  was  sufficient  to  cause  displacements  that 
could  be  easily  imaged  with  the  CCD  camera. 
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Figure  3.  (A)  Force  on  the  glass  fibre  in  N  versus  displacement  in  m.  (B)  Displacement  of  the 
Teflon  bead  attached  to  the  fibre  versus  average  electrical  power  delivered  to  the  transducer. 


A  3  mm  diameter  Teflon  bead  was  attached  to  a  glass  fibre  and  suspended  at  the  centre 
of  the  sample  cell,  which  was  filled  with  water.  The  glass  fibre  was  calibrated  by  attaching  it 
to  a  precision  translation  stage  and  using  a  laboratory  balance  to  determine  force  versus  fibre 
displacement.  The  data  shown  in  figure  3(A)  gave  a  force  constant  of  2.18  x  10“2  N  m-1  for 
the  fibre.  The  total  force  on  a  bead  suspended  from  the  fibre  is  the  sum  of  the  restoring  force 
from  the  fibre  and  the  force  from  gravity,  which  can  be  written  for  small  angles  of  deviation 
from  the  vertical  as 


exp 


glassfibre  ' 


^bead 


(8) 


where  rabead  is  the  mass  of  the  bead,  g  is  the  gravitational  constant,  L  is  the  length  of  the 
fibre  and  d  is  its  displacement.  For  the  experiments  described  here,  the  mass  of  the  bead  was 
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Figure  4.  (A)  Acoustic  pressure  in  Pa  from  the  calibrated  transducer  versus  driving  voltage  from 
the  power  amplifier.  (B)  Force  on  the  bead  in  N  versus  the  square  of  the  pressure  from  the 
transducer. 

35  mg,  and  the  length  of  the  glass  fibre  was  10  cm.  The  displacement  of  the  bead  versus  the 
time  averaged  electrical  power  delivered  to  the  transducer  is  shown  in  figure  3(B).  The  break 
in  the  slope  of  the  curve  can  be  attributed  to  a  departure  of  the  efficiency  of  the  transducer 
from  linearity  at  high  driving  voltages.  The  acoustic  pressure  generated  by  the  transducer 
at  various  power  levels  from  the  power  amplifier  was  determined  using  a  calibrated  LiNiC>3 
transducer  (LaserSonics  Technologies,  Model  WAT-04).  The  plot  of  acoustic  pressure  from 
the  calibrated  transducer  versus  average  driving  voltage  applied  to  the  transducer  is  shown  in 
figure  4(A)  which  gives  a  curve  with  the  same  qualitative  features  as  the  curve  in  figure  3(B), 
and  permits  calibration  of  the  transducer  pressure  in  terms  of  the  power  amplifier  output 
driving  voltage.  With  this  calibration,  the  plot  in  figure  3(B)  can  be  converted  to  a  plot  of 
force  versus  the  square  of  the  acoustic  pressure  as  shown  in  figure  4(B).  The  theoretical  curve 
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from  equation  (6),  also  plotted  in  figure  4(B),  shows  reasonably  good  agreement  with  the 
experimental  measurements.  The  maximum  force  recorded  in  the  experiments  was  5  /xN. 


4.  Acoustically  modulated  x-ray  phase  contrast  images 


The  method  of  modifying  a  phase  contrast  image  using  acoustic  radiation  pressure  (Bailat  et  al 
2004)  consists  of  two  steps:  first,  an  x-ray  image  is  made  with  a  sound  beam  directed  into  a 
body  to  displace  an  object  through  acoustic  radiation  force  and  the  image  stored  in  a  computer, 
second,  another  x-ray  image  of  the  object  is  taken,  this  time  without  the  presence  of  the  sound 
beam  and  the  image  is  recorded  in  the  computer.  The  two  images  are  then  subtracted  pixel  by 
pixel  to  give  a  subtracted  phase  contrast  image,  the  component  of  the  image  from  absorption 
contrast  being  largely  eliminated,  leaving  a  nearly  pure  phase  contrast  image,  inherently 
background  and  flatfield  corrected. 

The  principle  of  the  method  was  demonstrated  by  taking  images  of  two  3  mm  diameter 
Teflon  beads  cast  in  a  6.5  cm  long  block  of  agarose.  As  shown  in  figure  5(A),  the  two  beads 
were  both  within  the  field  of  the  x-ray  beam,  but  only  one  bead  was  irradiated  with  ultrasound. 
Figure  5(B)  shows  an  image  of  the  two  beads  where  the  microfocus  x-ray  tube  was  operated 
at  high  power,  80  W,  to  enlarge  the  diameter  of  the  x-ray  source  to  approximately  75  /xm. 
With  this  source  diameter,  the  phase  contrast  component  of  the  image  disappears  as  a  result 
of  convolution  of  the  phase  contrast  fringes  over  the  source  area  of  the  x-ray  beam.  The 
exposure  time  for  the  image  was  75  s.  In  figure  5(C)  a  conventional  phase  contrast  image 
of  the  two  beads  is  shown  using  the  same  total  x-ray  fluence  as  was  used  in  figure  5(B),  but 
with  the  x-ray  tube  operating  at  10  W,  where  the  source  diameter  was  25  /xm.  The  fringes 
at  the  perimeters  of  the  beads,  which  serve  to  define  the  edges  of  the  beads,  are  the  result  of 
interference  from  rapid  phase  variations  in  the  x-ray  paths,  which,  for  a  sphere,  are  largest  at 
its  perimeter  (Wilkins  et  al  1996).  The  results  of  the  image  subtraction  with  the  ultrasound 
directed  onto  only  the  bead  closest  to  the  x-ray  tube  are  shown  in  figure  5(D).  The  bead 
that  was  not  irradiated  with  ultrasound  is  not  visible  in  the  image,  while  the  phase  contrast 
component  of  the  image  of  the  irradiated  bead  shows  up  with  the  phase  contrast  component 
highlighted.  The  subtracted  phase  contrast  image  in  figure  5(D)  shows  a  reversal  of  the  bright 
and  dark  regions  at  the  perimeter  of  the  bead  in  the  direction  of  the  ultrasound  at  the  two  sides 
of  the  bead  (which  appear  at  the  top  and  bottom  of  figure  5(D)):  the  interior  of  the  arc  at  the 
top  of  the  figure  appears  bright  while  it  is  the  exterior  of  the  arc  at  the  bottom  of  the  figure 
that  appears  bright.  The  reversal  of  the  shading  is  a  consequence  of  subtraction — the  opposite 
shading  could  be  produced  by  reversal  of  the  order  of  subtraction  of  the  two  images. 

For  our  purposes,  an  analysis  of  the  ultrasound-induced  temperature  increase  in  the 
selected  tissue  is  important.  A  first  approximation  of  the  temperature  change  is  given  by 
2  /xv 

AT  =  —It 

y 

where  /x  is  the  ultrasound  attenuation  coefficient  in  soft  tissue  (=0.05  cm-1  MHz-1),  v  is 
the  ultrasound  centre  frequency  (=  1  MHz),  y  is  the  heat  capacity  per  unit  volume  of  soft 
tissue  (=4  W  s  cm-3  °C_1),  I  is  the  ultrasound  intensity  in  W  cm-2  and  t  is  the  application 
time  (15  /xs)  (Herman  and  Harris  2002).  The  maximum  intensity  can  be  calculated  from  the 
square  of  our  maximum  recorded  pressure  p2ax  using 

p 2 

/  =  (10) 
2  pc 


(9) 


where  /?^ax  =  1.5  x  1012  Pa2  (Morse  1981).  We  obtain  I  =  5  x  105  W  m  2  =  50  W  cm  2 
and  a  subsequent  temperature  increase  of  approximately  1.8  x  10-5  °C.  However  in  our 
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Figure  5.  (A)  Side  view  of  the  sample  consisting  of  two  Teflon  beads  cast  in  agarose.  The  two 
beads  are  separated  by  3  cm.  (B)  Image  of  the  beads  with  no  phase  contrast.  (C)  Image  of  the  beads 
with  phase  contrast.  (D)  Subtracted  image  of  the  two  beads.  The  bead  irradiated  with  ultrasound 
is  nearest  the  CCD  camera  and  produces  a  slightly  larger  image  than  the  bead  located  outside  the 
sound  beam.  The  x-ray  tube  voltage  for  the  images  was  90  kV.  The  image  in  (B)  was  taken  with 
a  75  s  exposure  time  with  the  x-ray  tube  operating  at  80  W;  the  image  in  (C)  was  taken  with  eight 
times  lower  tube  wattage  and  eight  times  longer  exposure  to  give  the  same  total  x-ray  fluence. 


case,  burst  ultrasound,  a  more  specific  model  which  considers  an  ultrasound-bone  interaction 
should  be  used  as  a  worst-case  scenario.  The  upper  temperature  limit  for  burst  ultrasound  is 
approximated  by 

AT  =  189.74  x  10"6/V7  (11) 

with  t  in  s  and  I  in  W  m-2.  Using  equation  (11)  with  the  same  data  given  above  results  in 
an  estimated  temperature  increase  of  0.37  °C.  It  should  be  noted  that  since  (11)  was  derived 
specifically  for  the  case  of  an  ultrasound-bone  interface  the  resulting  temperature  increase 
in  the  case  of  soft  tissue  would  actually  be  smaller  and  would  not  present  any  dangers  in 
biological  applications. 

In  order  to  evaluate  the  method  using  tissue-like  samples,  a  phantom  of  micron-sized 
agarose  objects  embedded  in  chicken  breast  was  examined.  The  objects  in  the  phantom 
consisted  of  pieces  of  a  casting  made  by  mixing  carbon  nanoparticles  with  agarose  to  achieve 
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Figure  6.  X-ray  images  of  carbon  enriched  agarose  objects  embedded  in  chicken  flesh:  the  round 
objects  are  fat  bubbles  in  the  chicken  flesh.  The  entire  sample  is  immersed  in  water.  (A)  and  (B) 
conventional  phase  contrast  images  taken  with  a  25  /xm  diameter  x-ray  source,  and  (C)  and  (D) 
corresponding  ultrasound  enhanced  x-ray  phase  contrast  image. 

a  density  approximately  14%  higher  than  that  of  the  chicken  breast,  comparable  to  the  density 
difference  between  malignant  and  healthy  human  breast  tissue  as  reported  in  Duck  (1990), 
Johns  and  Yaffe  (1987).  The  conventional  phase  contrast  images  in  (A)  and  (B)  for  two 
different  regions  of  the  phantom  show  fat  globules  as  being  clearly  outlined  elliptical  objects 
and  the  agarose  objects  as  the  dark  masses,  their  estimated  sizes  being  of  the  order  of  30  x 
30  /xm  in  (A)  and  30  x  60  /xm  in  (B).  The  ultrasound  contrast-enhanced  x-ray  images 
corresponding  to  (A)  and  (B)  are  shown  in  (C)  and  (D),  respectively.  Both  regions  of  the 
phantom  were  irradiated  with  ultrasound.  The  subtracted  images  are  seen  to  have  increased 
contrast  of  the  agarose  objects,  and  the  fat  globules  as  well,  highlighting  those  parts  of  the 
phantom  displaced  by  the  ultrasound.  Both  the  fat  globules  and  the  higher  density  agarose 
objects  are  surrounded  with  identical  shading,  indicating  motion  in  the  same  direction  under 
the  influence  of  the  acoustic  radiation  pressure.  Some  of  the  objects,  well  defined  on  the 
images  (A)  and  (B),  are  no  longer  visible  in  subtracted  images  and  can  therefore  be  identified 
as  either  artefacts  arising  from  optical  imperfections  in  the  sample  container  or  imaging  optics. 
Note  that  a  simple  differential  shift  of  two  phase  contrast  images  followed  by  their  subtraction, 
as  opposed  to  the  present  acoustically  modified  procedure,  would  leave  artefacts  from  optical 
imperfections  in  the  final  image. 


5.  Discussion 

The  intensity  profile  of  the  subtracted  image  S(x)  can  be  written  as  S(x)  =  I  (x  +  h)  —  I  (x) 
where  h  is  the  displacement  caused  by  the  application  of  ultrasound.  For  small  h ,  it  is 
appropriate  to  identify  the  subtracted  image  as  approximately  proportional  to  the  space 
derivative  of  the  phase  contrast  image  in  the  direction  of  propagation  of  the  ultrasound, 
^(jt)  =  hf(x).  As  such,  low  spatial  frequency  features  such  as  broad  absorptions  are 
suppressed  in  the  image  while  the  high  frequency  features  such  as  those  at  interfaces  between 
dissimilar  materials  are  emphasized.  When  the  acoustic  displacement  of  an  object  is  of 
the  order  of  the  fringe  spacing  of  the  phase  contrast,  the  highlighting  of  the  phase  contrast 
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component  of  the  image  is  maximal;  the  image  of  the  object  is  then  almost  a  pure  phase  contrast 
image  responsive  to  changes  in  density.  In  general,  the  contrast  in  the  image  is  determined 
by  the  combined  effects  of  the  tissue  selectivity  of  the  ultrasound  based  on  differences  in 
acoustic  impedance,  together  with  the  response  of  phase  contrast  to  density  variation.  The 
resolution,  on  the  other  hand,  is  governed  by  the  parameters  of  conventional  phase  contrast 
imaging,  which  include  the  degree  of  spatial  coherence  of  the  x-ray  beam  and  the  sample-to- 
source  and  sample-to-image  distances.  A  distinctive  feature  of  the  method  described  here  is  its 
vector  character — the  subtracted  image  shows  only  those  features  of  the  object  moved  along 
the  propagation  direction  of  the  ultrasound,  which  is  evident  in  figures  5  and  6.  Clearly,  it  is 
possible  to  make  multiple  images  where  sound  is  directed  from  different  angles  into  a  body  and 
to  process  the  images  to  record  displacement  in  more  than  one  direction.  Additionally,  there 
are  different  ultrasound  excitation  schemes  that  can  be  employed  leading  to,  for  instance, 
standing  acoustic  waves,  that  can  be  generated  within  a  body  to  give  contrast  dependent 
on  eigenmodes  of  vibration  and  other  elastic  properties  of  the  body,  further  enhancing  the 
capability  of  x-rays  to  probe  subtle  physical  properties  of  matter. 
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